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1.0 SUMMARY 


Program L225 (TEV126) calculates the time responses of a linear system by convoluting 
the impulsive response function with the time-dependent excitation. Fast Fourier 
transform (FFT) techniques are used to perform the convolution in the frequency 
domain and transform the results into the time domain. 

The program accepts the frequency response functions of a linear system, multiplies 
them by the Fourier transform of the desired time-dependent excitation function, and 
performs the inverse Fourier transform of the product to determine the response time 
histories. The Cooley-Tukey algorithm is used to evaluate the Fourier transforms. The 
analysis procedures are subject to the following theoretical limitations: 

• System being considered is linear. 

• System is stable. 

• System is initially in steady state equilibrium. 

• Time response of the system is real. 

• Time response is negligible beyond the maximum time of interest. 

• Forcing function begins at time zero, is real and is negligible beyond the maximum 
time of interest. 

Data are input to the program on both cards and magnetic files (tapes or disks). The 
user has the capability of using one of seven standard excitation functions contained in 
the program or of defining an arbitrary excitation function on cards or magnetic file. 

Output from L225 consists of: 

• Printed time histories 

• Maximum and minimum values of the time responses 

• Plotting vectors of the time histories 

• Plotting vectors of convoluted responses in the frequency domain 


The significant program limitations that restrict the size of the problem being analyzed 
are: 

• The maximum number of response functions, either generalized coordinate 
responses or load responses, is 100. 

• The maximum number of frequencies at which the input response functions are 
defined is 250. 

• All interpolation of frequency response functions is linear. 


2.0 INTRODUCTION 


The computer program L225 (TEV126) was developed for use as either a standalone 
program or as a module of a program system called DYLOFLEX developed for NASA 
under contract NAS1-13918 (ref. 1). Because of the DYLOFLEX contract requirements 
defined in reference 2, a program was needed to calculate the time histories of different 
airplane responses to discrete excitations using frequency response data calculated in a 
power spectral density (PSD) solution. An existing program 1 ’ was revised according to 
established specifications’* to make it compatible with other modules of the DYLOFLEX 
system. 

The purpose of this volume is to aid the person wishing to use the program to calculate 
response time histories. To meet this objective, the following items are discussed in 
detail: 

• Theoretical formulation of the problem 

• Design and structure of the program 

• Guidelines to the actual use of the program 

A sample problem is also presented in this volume to aid the user in execution of the 
program. 


^mith, B. E.; and Johnson, M. R.: Discrete Gust Convolution Program (TEV126). Boeing 
document D6-29668TN, January 1969. 

2 Every, J. S.: Frequency Response Function Sorting Program (TEV126, LINK1, DIG1). 
Boeing document D6-29669TN, vol. II, April 1971. 

3 Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System 
to Interface With FLEXSTAB. NASA contract NAS1-13918), BCS-G0701, September 1975. 
(Internal document) 


2 


I 


3.0 SYMBOLS AND ABBREVIATIONS 


The following 
document. 

Engineering 

notation 

f 

F a 

FFT 

f HARM 

f MAX0 

f maxi 

f(t) 

F(to) 

h(t) 

H(tu) 


list contains items that appear in the discussion sections of this 
Definition 

Cyclic frequency, cps. 

Aliasing frequency, cps. 

Fast Fourier transform 

Maximum frequency for which the FFT is defined, cps. 

Maximum frequency for which the input frequency response functions 
are considered defined, cps. 

Cutoff frequency, cps. 

Forcing function expressed in the time domain. 

Forcing function expressed in n, the frequency domain. 

Impulsive response function. 

Input frequency response function (complex). 


M 

N 

t 

t max 

x(t) 

X(«) 

y(t) 

Y(<u) 

Y *(to) 


Power of 2 which is used to define the number of points to be used in 
the FFT. 

Number of points used in the FFT. 
time, sec. 

Maximum time for which the response time history is calculated, sec. 
Response quantity expressed in the time domain. 

Response quantity expressed in the frequency domain. 

General time function. 

General frequency response function (complex). 

Complex conjugate of a general frequency response. 


3 


Af 


Frequency interval, cps. 


At 

CD 

T[f(t)] 


Time interval, sec. 

Angular frequency, rad/sec. 

Fourier transform of the forcing function, fTt). 



4.0 ENGINEERING AND MATHEMATICAL DESCRIPTION 


4.1 MATHEMATICAL FORMULATION 

The motion of a linear system can be defined by a set of linear differential equations. 
The response of the system to a particular type of force is obtained from the solution of 
these equations with respect to that force. Of particular interest are the responses of a 
linear system to two particular types of forces: a harmonic excitation and an impulsive 
excitation. 


The steady state response of the system (the response of the system after all transients 
have died out) to harmonic excitation is described by the function H(a>), called the 
frequency response function. This function is used extensively in airplane power 
spectral analysis and in dealing with frequency domain aerodynamics. 

The solution of the system’s differential equation to a unit impulse is described by the 
impulsive response function h(t). The value in the impulsive response function lies in its 
use with the principle of superposition, which is only valid for linear systems. Any 
arbitrary forcing function can be considered as an infinite set of impulses of varying 
magnitude. The responses of the system to the impulses can be added together to give 
the total response due to the arbitrary forcing function. This superposition process is 
mathematically written as the convolution integral 


x(t) = / 


t 

o 


h (t - r) f (r)dr 


(1) 


where f(t) is the arbitrary forcing function (ref. 3). 

Therefore, having H(a>) and h(t) defined for a system, the response to both continuous 
and discrete excitation is defined. It would be advantageous, however, to be able to use 
one of the functions, H(cu) or h(t), to determine the other. Under certain conditions, 
Fourier transform analysis accomplishes this task. 


A function, y(t), defined in the time domain, can be transformed or completely described 
in the frequency domain by applying the transform 


Y(a») = / °° y(t)e' lajt dt 
-00 


(2) 


Conversely, if Y(w) is defined, the inverse transform can describe the function in the 
time domain. 


Y(t) 


-tJ' 


Y(&>)e lb>t dcu 


(3) 
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These integrals exist provided the following condition is met, that 


J 


y(t) 


dt 


is finite. If this is so, then y(t) and Y (co) are known as Fourier transform pairs. 


( 4 ) 


For a stable system, Fourier transform theory (ref. 3) shows that the impulsive response 
function h(t) and the frequency response function H(a>) are transform pairs; i.e., 


H(ftj) = / °° h(t)e' i£ot dt (5) 

-00 

and 

h(t) = — f H(w) e lfl>l dcu (6) 

277 -00 

Therefore, once the system’s response is determined by either h(t) or H(o>), the other 
function can be calculated. 

Besides providing a link between the time domain and the frequency domain for a 
stable linear system, Fourier transform theory facilitates the evaluation of the 
convolution integral given by equation (1). In the frequency domain, convolution 
becomes a simple multiplication of functions (ref. 4), so that 

X(«u) = H(w) F(oi) (7) 

where X(co), H(o>), and F(o>) are transform pairs of x(t), h(t), and fit), respectively. 

Using the transform relationships and the frequency domain representation of the 
convolution integral, it becomes possible to develop the response of a linear system in 
the frequency domain and use the results to determine system’s response to a discrete 
time function. 

4.2 DIGITAL APPROXIMATION OF THE FOURIER INTEGRAL 

The mathematical formulation dealt with continuous functions and infinite limits of 
integration. Digital computers are incapable of dealing with either, therefore; 
approximations must be introduced to be able to deal with Fourier transform integrals. 

First, by placing restrictions on the type of functions which can be dealt with, the effect 
of introducing finite limits into the integrations of equations (2) and (3) can be 
minimized. The restrictions imposed on y(t) are; 

• y(t) must be real 

• y(t) must be zero for t < 0 

• y(t) must be negligible for t > t MAX 



The last two restrictions enable equation (2) to be written with finite limits. 

Y(tu) = / TMAX y(t) e‘ ift>t dt (8) 

0 

The first restriction, y(t) must be real, results in the frequency domain property of 

Y(-w) = Y *(o>) (9) 

where Y *(to) is the complex conjugate of Y(to). This property coupled with the 
restriction that Y(<o) is negligible for ui > 27rFjy{AXl enables the approximation of the 
inverse relationship by the equation 


y(t) = 2Real { f 2 ” FmaX1 Y(a))e' la,t dt } (10) 

0 

The relationship between the maximum time T\jax an< i the maximum frequency 
FmaXI w >11 be discussed in section 4.5. 

The second variation from the exact formulation, which must be considered when 
dealing with digital computation, is the use of digitized functions rather than 
continuous functions. Using a function described by a discrete number of points is 
mathematically equivalent to time sampling a continuous function. Time sampling 
introduces an aberration into the transform process which muBt be taken into account. 
Figure 1 illustrates a time function y(t) and the Fourier transform Y(f). 




(A) (B) 


Figure 1.— Function in Time and Frequency Domain 


The sampling of y(t) may be regarded as multiplying y(t) by a series of unit impulses at 
intervals At. Figure 2 shows a series of impulses and the Fourier transform pairs. 
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Figure 2.— Impulse Functions in Time and Frequency Domain 

From Fourier transform theory, the transform of the product of two functions is the 
convolution of the transform pairs, or 

co -ia>t 

f y(t)yi(t)e dt = Y(o))*Yi(o>) (11) 

-OO 

The graphic result of equation (11) is shown in figure 3. 



The result of sampling then is to make the transform function periodic with a period 
equal to 1/At. This results in an overlapping or aliasing at the frequency F a = l/2At; 
therefore, caution should be exercised in choosing the size of the sampling interval At. 
Sampling in the frequency domain results in aliasing in the time domain. A further 
discussion of aliasing appears in section 4.2.2. 

The digital forms of equations (8) and (10) are 


N-l 

Y(jAf) = At ^I^y(kAt) e 


- i 2ir (jAf kAt) 


where j = 0, 1, . . ., N-l and 


i i27r(jAf kAt) I 

Y(kAt) = 2Real Af L Y(jAf) e 

I j = 0 ' 


( 12 ) 


(13) 
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where K = 0, 1, . . N-l and N is the number of sample points. Equations (12) and (13) 
were derived using the trapezoidal rule approximation. 

The condition that the time response is zero for t<0 allows equation (13) to be evaluated 
by two independent methods. These can then be compared as a check on the correctness 
of the solution. Rewriting equation (13) 


N-l 

Y(kAt) = 2Af E [YftQAf) cos(2irjAf kAt) - iYjfjAf) sin(27rjAf kAt)] (14) 

j=0 

where y(t) can be seen to be the sum of an even and an odd function of time. These add 
for t>0 and for t<0; they subtract to yield a zero result. Therefore, y(t) can be evaluated 
from either the real or the imaginary part of Y(f); i.e., 

N-l 

Y(kAt) = 4Af E YR(jAf) cos(2irjAfkAt) (15) 

j = 0 


N-l 

Y(kAt) = 4Af E Yi(jAf) sin(27rjAfkAt) (16) 

j = 0 

4.3 ANALYSIS PROCEDURES 

The purpose of the L225 (TEV126) program is to evaluate the convolution integral 


x(t) = / T MAX 
3 0 


h(r)f(t-r)dT 


(17) 


using Fourier transform techniques. In the frequency domain, equation (17) is written 
as 

X(o>) = H(qj)F(o») (18) 


The time history program accepts both the frequency response function and the forcing 
function defined in the time domain and performs the necessary transformations, 
multiplications, and inverse transformations. 


Certain restrictions are assumed with regard to the impulsive response function h(t) 
and the forcing function f(t): 


h(t) = 


f(t) = 


0 

Nonzero 

Negligible 

0 

Nonzero 

Negligible 


t<0 

0=£t«T M AX 

T <t 

J MAX 

t<0 

O^t^TjyjAX 

T MAX <t 
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The restrictions are a result of the desire to minimize any error introduced by the 
approximations made in evaluating the transform integrals. (This last restriction 
concerning h(t) implies h(t) corresponds to a stable system.) 

The process of obtaining the time response begins with the input of the frequency 
response function H(o). The program can accept frequency response functions derived in 
either one of two ways. First, if the solution to the differential equation was of the form 
e la> , then this is equivalent to the impulsive response function - frequency response 
function relationship 


. 00 -iot 

H(o) = / h(t) e dt 
-00 


(19) 


Second, if the solution assumed was of the form e ia)t ’ then the equivalent h(t) - H(o) 
relationship is 


. icut 

H(cu) = J h(t) e dt (20) 

-00 

The program assumes the form given by equation (19); , therefore, if the e lwt form of 
the solution was used to derive H(a>), the user can exercise the option to use the complex 
conjugate of H(a>) for the time history analysis. Thus, H(a») will be defined consistently 
with the program assumption. The maximum frequency for which H(o) is defined is 
denoted as F^axO- ^ i® assumed that the value of H(oi) beyond Fmaxo i® negligible. 
The program also assumes that the real part of H(o) is an even function and that the 
imaginary part of H(co) is an odd function. 

Input of H(o) may be at values of other than those used in the transform equations (15) 
and (16). Linear interpolation is used between input frequencies to obtain the values of 
H(o) at the required integration frequencies. It is important to input H(a>) with a 
sufficient density to allow for accurate linear interpolation. 

After defining the frequency response function, the next step is to determine the 
Fourier transform of the forcing function. The program has a number of internally 
defined forcing functions and also haB the capability of accepting arbitrary functions via 
external input. With fit) defined, the transform equation 

Tmay -icot 

F(w) = f f(t)e dt (21) 

0 


is numerically evaluated. 

The multiplication of equation (21) is performed, and the inverse transform 

2*-Fm A xi 


x 0 (t) = 2Real { — — / 


x(o)e do | 


(22) 


10 



is evaluated. Because of the restrictions placed on h(t) and f(t), the time response should 
be zero for t<0 and a real function. To check the accuracy of the numerical integrations 
and the validity of the assumptions, two alternate expressions may be used to evaluate 
x(t) 


4 . 2 ttFmAX 1 

xj(t) = J Real{x(co)} coso>tda> (23) 

2tt o 


4 27rF MAXl . . 

X 2 (t) = - — j Imaginary{x(co)} sin<otda> (24) 

27T o ' 

The assumption that the contribution of X(o>) to x(t) beyond the cutoff frequency F^^xi 
is negligible has also been included in the expressions. Notice that xo(t) is the average 
of x^(t) and X 2 (t) which would all be equal if the program error was zero and the 
physical restrictions imposed on h(t) and f(t) were adhered to exactly by the program 
input functions. 

Equations (21), (18), and (22) respectively, appear in digital form as: 


N-l 

T [ f(j At)] k = At f(jAt)e 


-ikAtojAt 


-At 


f(0) 


+ f((MM-l)AtV 


-ikAo>(MM-l)At 


X(kA<o) = H(kAo>)T [f(jAt)] k 
( i N-! 

xo(jAt) = 2 Real {- — Ao> Yi X(kAo>)e 

l 277 k=0 


ikAojjAt 


.i-Ao m 

2tt \ 2 


+ X((MM-l)A<u), 


>)f 


ijAt(MM-l)Aoj) 


(25) 

(26) 


(27) 


where: 



j 

= 0, 1, . . . 

, N-l 

k 

II 

O 

H-* 

, N-l 


( Nonzero 

k = 0, 1 NN 

x(kAa>) 

= 



V Zero 

k = NN + l, . . ., N-l 
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f(JAt) 

= - 

I’ Nonzero 
[ Zero 

j = 0, 1, . . ., NN 
j = NN + N-l 

MM 

= 

NN+1 

if NN=sN-l 

MM 

= 

NN if 

NN=N 


The evaluation of equations (25), (26), and (27) is made using both the Cooley-Tukey 
and fast Fourier transform algorithms (ref. 4). 

4.4 PROGRAM ACCURACY 

The program accuracy is a function of the program input functions, the cutoff frequency, 
and the spacing in the time and frequency domains. For optimum conditions there 
should be no discontinuities in the input functions and appropriate spacing should be 
used in the time and frequency domains. However, the use of the time history for a 
variety of problems necessitates the user’s understanding of the effects of the FFT 
approximations. Reference 4 gives a thorough explanation of FFT techniques and the 
impact of the various approximations made in the FFT process. A brief discussion of 
how to set up a time history solution follows. 

4.4.1 CUTOFF FREQUENCY CONSIDERATIONS AND EFFECTS 

In the program analysis, the product of the frequency response function H(f) and the 
transform of the forcing function T[f(t)] is assumed to be symmetric about the origin 
and negligible for frequencies beyond the cutoff frequency F\ia.X 1 ( 8ee fig- 4). 

X (f ) = H(f) T[f(t)] 



Frequency 

F igure 4. — Transform of the Time Response 

The effects of truncation in the frequency domain are a superposition of waviness on the 
original time function x(t> and a possibility of causing small values of x(t) at t = 0. The 
amount of waviness and the magnitude of the initial value will depend on the type of 
function being considered and the value of F^AXl- Functions with low frequency 
content, such as sine and cosine functions and gradual ramps, will show little effect of 
truncation. However, functions having high frequency content, such as steps and sharp 
ramps, will be affected more by the truncation function. Thus, it is important to pick a 
value of Fmaxi that will include as much of the frequency function as possible. 


12 


When considering a frequency response function and/or a particular excitation function 
for the first time, the user should exercise the option to obtain a plot of | X(f) | , where 
|X(f)| is the magnitude of the complex function X(f). The range of f considered should 
be as high as possible while keeping in mind the frequency response function and 
forcing function used. This plot will give an indication as to what value of cutoff 
frequency will give suitable results. If the algorithm is to give accurate results, the 
values of |X(f)| beyond the cutoff frequency should be negligible. 

It is important to note that the approximation being introduced by assigning a cutoff 
frequency is that 


/ co i2-n-ft . 

E(t,F) = 2Real { J X(f)e df 
( F I 


is negligible when compared to 


1(T,F) = 2Real { J F X(f)e l27rft df } 

• 0 ) 


(28) 


(29) 


that is 


x(t) = I(t,F) + E(t,F) 
= Kt,F) 


(30) 


Looking at the plot of | X(f) | does not give an absolute indication as to whether or not a 
particular value of cutoff frequency will give satisfactory results. It should only serve as 
a guide. 


To illustrate the effect of cutoff frequency, an input frequency response function of a 
constant value of one was convoluted with two forcing functions: a rectangular wave 
and a one-minus cosine wave. Three values of cutoff frequencies were considered, 2.5, 
10, and 20 cps. Figure 5 shows the rectangular wave and how it is affected by the cutoff 
frequency. Figure 6 shows the effect of cutoff frequency on a one-minus cosine function. 
The two figures demonstrate that, for functions with sharp discontinuities, high values 
of cutoff frequency are needed, whereas, smooth functions do not require as large a 
frequency range. Functions with sharp discontinuities have high frequency components, 
whereas, smooth functions contain little contribution from frequencies beyond the 
fundamental frequency of the wave. This can be seen from the plot of the frequency 
functions for the two waves shown in figure 7. The rectangular wave contains 
significant contributions from frequencies above 10 cps. The one-minus cosine wave, 
however, shows little or no contribution beyond 10 cps. Thus, the accuracy of the 
calculated time response of a function will depend on its frequency content and whether 
or not the frequency range considered adequately accounts for the frequency content. 
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TIME (SEC) 


Effect on a 




















Figure 7.— Frequency Response Functions for a Rectangular 
and a One Half ( 1 - cos) Wave 
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4.4.2 TIME AND FREQUENCY INTERVAL SIZE 


Accuracy of the results obtained is also a function of the time and frequency interval 
size chosen. The interval size is dependent on the number of points, N, chosen to be 
used in the FFT, the maximum time of interest, T\iax> and the maximum frequency 
which the input frequency response functions are considered defined, F^jaXO- The 
relationships of the variables will be discussed in section 4.5. It is the purpose of this 
section to illustrate the effect that interval size has on the results obtained. 

It was mentioned in section 4.2 that digitizing a continuous function in one domain 
(time or frequency) results in an aliasing (or overlapping) effect in the other domain. 
This aliasing is directly related to interval size. As shown in figure 3, the sampling of a 
time function results in its transform point being repeated in the frequency domain at a 
frequency 


F a = 1/At (31 ) 

Similar sampling of a frequency function will result in a time repetition at 

T a = 1/Af (32) 

To minimize the distortion of the frequency and/or time functions, the interval sizes 
should be such that the portion of the frequency function beyond F a /2 or the portion of 
the time function beyond T a should be negligible. Note that for the frequency functions, 
half F a is the point of overlap because the frequency functions are symmetric about the 
origin (see fig. 3). On the other hand, the time functions have no values for t < 0, 
therefore, T a is the point of overlapping. 

To illustrate the aliasing effect, time responses of the first elastic mode and wingtip 
acceleration to a one-minus cosine forcing function were calculated using two different 
sample sizes. The first case used a frequency interval size of 0.0395 which gives an 
aliasing time T a = 25.3 sec. The second case used an interval size of 0.289 which gives 
an aliasing time of 3.5 sec. Figures 8 and 9 show the time histories of these response 
quantities. Case one results show that the response of both quantities damp out beyond 
8 sec, therefore, overlapping of the time function beyond 25.3 sec has no effect on the 
time history. However, both quantities show significant response between 3 5 and 
7.0 sec. It would be expected, therefore, that any overlapping that occurs within this 
time frame would significantly distort the resulting time history. The results of case 
two show this to be true. 


The analysis procedures used by L225 require the analyst to account for aliasing effects 
in both the time and frequency domain. Aliasing in the frequency domain must be 
considered with regard to the forcing function, f(t), since it is this function that is 
transformed from the time domain to the frequency domain. The choice of the size of the 
time sampling interval, At, must be made with regard to the nature of the forcing 
function. For example, a forcing function with sharp discontinuities may have 
significant high frequency content in its Fourier transform and thus may require a 
small time interval to ensure a large aliasing frequency. The frequency sampling size, 
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4 5 6 

TIME (SEC) 


TIME (SEC) 


Figure 8.— Interval Spacing Effect on First Elastic Mode Response 







Af, on the other hand, should be chosen with regard to the response quantity, x(f). High 
values of x(f) in the low frequency range may indicate that the resulting time function, 
x(t), could have significant magnitude at large values of time. Therefore, the size of the 
frequency interval should be such that the aliasing time, T a , is well beyond the 
maximum time of interest, T^aX- 

4.4.3 SAMPLE RESULTS 

As an illustration of the program accuracy, the results of a transport airplane analysis 
are presented. The airplane was subjected to three different forcing functions: a 
25-chord one-half (1-cos) function, a 25-chord rectangular wave, and a gust digitized 
from record flight data. Shear at one of the wing stations was chosen as the response 
quantity to be presented. Figures 10 through 12 show the magnitude |X(f)| of the 
frequency function of the convoluted shear. Figures 13 through 15 show the resulting 
time histories. Two different values of cutoff frequency were used. Note that the effect 
of cutoff frequency is more pronounced in the cases of the rectangular wave and the 
measured gust because in both of these cases there is significant response above 3.5 cps. 
Included in the time history figures are the results obtained using an initial value 
routine (a Runge-Kutta numerical integration method). 

4.5 ANALYSIS PARAMETER SELECTION 

As mentioned in section 4.4.2, the accuracy of results of an analysis will depend on 
the type of input functions (both forcing and response), the frequency and time interval 
spacing, and the cutoff frequency. The capability exists with card set 9.0 for the user to 
select one of several options to be used to calculate the analysis control data that will 
best suit the problem at hand. 

The analysis parameters are calculated based on the following relationships 

N = 2 M FharM = N*Af T MAX = NN *At (33) 

F HARM = 1/At F MAX1 = NN *Af 


where: 

N 


M 


F HARM 


At 


Af 

F MAX1 

t max 

NN 


Number of points to be used in the Cooley-Tukey algorithm for the 
FFT calculation 

Power of two used to define N 

Largest frequency value for which the Fourier transform is defined 
Time interval 
Frequency interval 
Cutoff frequency 

Maximum time for which the time response is calculated 

The number of points needed to achieve integration to FmaXI or 

Tmax 
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Figure 10. —Airplane Response to a 25- Chord Rectangular Wave 
in the Frequency Domain 
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Figure 13.— Comparison of an Airplane Response to a 25-Chord 
Rectangular Wave in the Time Domain 
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Several options exist that enable the user to determine the parameters given in 
equation (33). All options require the user to define FmaxO- Table 1 is a list of the 
options available; the input data required with each option; and the restrictions under 
which each option operates. 


Table l.—List of Analysis Parameter Options 


Option 

Input data 

Restrictive conditions 

1 

None 

F HARM =F MAX1 =F MAX0 
M ( 1 ) = 10 

t max =1024/f maxo 

At=1/F MAX0 

2 

M 

F HARM =F MAX1 = F MAX0 

t max =n/f maxo 
A t= 1/f maxo 

3 

m, t 

f harm =f max ,<f maxo 


'max 

t max > n/f maxo 

At>1/F MAX0 

4 

M,At 

F HARM = F MAX1 <F MAX0 
t max >N/ f max0 

At>1/r 



r MAX0 



F HARM =F MAX1 =F MAX0 
IF T MAX < ' N ^ F MAX0* usual situation ) 

5 

M.At 

At<N/NN*F MAXQ 

or 

6 

m - t max 

F HARM =F MAX1 <F MAX0 

,F t max >n/f maxo 



At>N/NN*F MAXQ 

7 

M ' T MAX' At 

F HARM >F MAX1 



f maxo >f maxi 

At is changed only as a last 



resort 

8 

M ' T MAX' At 

F HARM >F MAX1 

f maxo >f maxi 

Tmax ' s changed only as a last 




resort 
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The parameters are defined under the two basic restrictions that 


F MAX1 55 F HARM (34) 

and 

f MAX 1 € f MAXO (35) 


where 

f MAXO = Maximum frequency for which the input frequency response functions 
are considered defined. 


Since the frequency function X(&j) is formed from the product of the input frequency 
response function, H(<u), and the Fourier transform of the forcing function, F(a>) (eq. 18), 
it makes no sense to integrate beyond the point where either F(a>) or H(oj) is adequately 
defined. The program restrictions stated in equations (34) and (35) result in such a check 
being present during parameter selection. If the user were to select a combination of 
parameters that would violate these restrictions, the program will automatically 
readjust the parameter so that the restrictions of (34) and (35) are met. The program 
does not, however, have any safeguards with regard to aliasing effects. It is the 
responsibility of the user to ensure that the selection of the analysis parameters given 
in table 1 will not introduce serious aliasing effects with regards to F(gj) or x(t). 


Options 1 through 4 are generally used to define the response for large values of time. 
When the maximum time of interest is small yet a large number of points are needed to 
properly define the transform, options 5 and 6 are recommended. These two options are 
normally selected in the majority of analyses. In options 7 and 8, the cutoff frequency is 
calculated from the relationship 


F MAX1 = T MAX /At2 * N 


(36) 


which is derived from equation (33). 

In options 3 through 6, the input variable At or T^AX * 8 subjected to change if the 
restrictive conditions are not met. The user has the option to have such changes printed. 
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5.0 PROGRAM DESIGN AND STRUCTURE 


The time history program L225 (TEV126) consists of three overlays. Figure 16 shows 
the overlay structure used for this program. The main overlay is a small program used 
to drive either or both of the two remaining overlays. The function of the SORT overlay 
is to select and rearrange the frequency response functions input from the PSD 
program. Input to this overlay is via cards and magnetic files. Output from SORT is the 
magnetic file "SRTFRP” containing the chosen frequency response functions. 


The ANAL overlay performs all convolutions and transformations. Card and magnetic 
file inputs are used. The user has the option to output printed time histories, a file of 
plotting vectors for the response time history, and/or a file of plotting vectors for the 
magnitudes of the output frequency response function | X(f) | . 
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Figure 16. — L225 (TEV126) Overlay System 
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6.0 COMPUTER PROGRAM USAGE 


6.1 MACHINE REQUIREMENTS 


The machine requirements for an execution of the program L225 (TEV126) are: 

Card reader Read control cards and data cards 

Printer Print calculated results 

Disk storage Contain a temporary file of sorted frequencies and 

responses 

Magnetic tape drive Read the program file, frequencies and responses, time 

functions, and write plot data files 


6.2 OPERATING SYSTEM 

Program L225 (TEV126) is written in FORTRAN for the CDC 6600. It may be compiled 
with either the RUN or FTN compilers. L225 (TEV126) may be executed under the 
KRONOS 2.1 or NOS operating system. 


6.3 CONTROL CARDS 

L225 (TEV126) may be extracted from a master file and executed on the KRONOS 2.1 
operating system with the following control cards and deck setup. 


Job Card 
Account Card 


REQUEST ( MASTER, F = I, LB = KL,VSN = 66XXXX) 
REWIND (MASTER) 

SKIPF (MASTER) / Retrieves the 

COPYBF (MASTER, L225) | absolute binary 

| overlays of L225 
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RETURN (MASTER) 


L225. 

EXIT. 


Prepare any input magnetic files 


Save any output magnetic files 


DMP(1,1 15000) 
— End-of-record 


l // Card Input 
Data 


— End-of-file 


6.4 RESOURCE ESTIMATES 
Field Length 

The field length required to load L225 (TEV126) is 100000 octal words. The minimum 
field length required to run the program is 115000 octal words. More may be required if 
the problem size is large. Estimate the field length needed with the following formula: 

FL(field length) = 20000 10 + 2(NFREQ + NRMAX) 


where: 

NFREQ = Number of frequencies 

NRMAX = Number of responses to be read; larger of NDOF and NLD (see card 

set 3.0) 

Note that FL must be changed to an octal base number before being placed on the job 
card. 


Timing 

The central processing time required to convolute one frequency response is 
approximately 2.5 sec. 

Printout 

The number of lines printed depends on the number of convolutions processed and the 
type of output requested. See cards 10.0 through 10.4. The maximum lines per 
convolution is approximately 2050. 
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Tape Driver 


A magnetic tape drive will be required if the program is to be accessed from a master 
tape, or if the frequencies and responses are on tape, or if the time functions are to be 
read from tape. 

Disk Storage 

The disk storage requirement of L225 (TEV126) is small enough that the user need not 
be concerned about it. The scratch file "SRTFRP” written in SORT, overlay (L225,l,0), 
will contain a maximum of 55,000 words. For the length of the plot files written by 
ANAL, overlay (L225,2,0), see section 6.6.2. 

6.5 INPUT DATA 

Program L225 (TEV126) reads input data from two sources: cards and magnetic files. 
Program directives and constants are input via cards. The frequency responses to be 
processed may be input on either cards or magnetic files. 

6.5.1 CARD INPUT DATA 

Figure 17 provides an overall picture of the card input data flow. The content and 
format of each card are described in detail in section 6.5.2. 

The task performed by L225 (TEV126) is broken into two subtasks, each with its own 
section of code known as a primary overlay. The two primary overlays are driven by a 
small program, the main overlay L225. 

L225 reads program directive cards to: 

• Assure that the data being read is intended for L225 (TEV126) 

• Determine which section of code (primary overlay) is to be executed next. 

• Terminate execution 

Each primary overlay has its own card input data sets. The primary overlays and data 
set numbers are listed in table 2. 

The following cards are all program directive cards read by L225. Note that they all 
begin with a keyword having "$” as the first character. 


Card 

Keyword 

Card 

Keyword 

1.0 

$TIME 

2.0 

$SORT 

1.1 

$TITLe 

7.0 

$ANALysis 

1.2 

$CHECk 

12.0 

$QUIT 


Each card set is described in section 6.5.2. 
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Card Set 7.0 
"$ ANAL" 



3e = NON- Stnd 



Card 8.4 
"TAPE " 




























Table 2. —Primary Overlays and Data Set Numbers 


Overlay 

Data 

Purpose of primary overlay 

Number 

Name 

sets 

1 

SORT 

3 through 6 

Sort frequency responses read from 
magnetic file or cards and write on 
a scratch file those that are to be 
retained for analysis in ANAL. 

2 

ANAL 

8 through 1 1 

Convolute the frequency responses 
read from the scratch file (prepared 
by SORT) and generate time histories. 


6.5.2 FORMAT OF CARD INPUT DATA 

All card data are read in fixed fields - specific columns of the cards. The required card 
columns are defined next to each keyword or variable on the pages that follow. Note the 
following conventions are used throughout the program: 

• All floating point variables are read with format E10.0 

• All integer variables are read with format 15. 

• All hollerith variables (keywords, etc.) are read with format A10 

Therefore, all data fields end on a card column which is a multiple of five. 

When the program is trying to recognize keywords, it checks only the first five 
characters. Any additional characters are ignored. 

The main overlay of L225 (TEV126) reads cards to determine which program module 
(SORT or ANAL) is to be called into execution. 

Note: Only the first five characters of any keyword are read to determine what is 

to be done. All underlined capital characters contained in the 
KEY WORD/VARI ABLE field of the input card sets must be left justified 
and punched in the card columns specified in the COLS, field of the input 
card sets. 
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Card Set 1.0-Program Directive 


' 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

ST I 'IE 

A10 

This card must always be the first card for the Time History 




L225(TEV126) program. 




Note: For any card with only a keyword on it, the user may use 




card columns 6-72 for his comments. 


Card 1.1 -Title Card 



KEYWORD/ 



COLS . 

VARIABLE 

FORMAT 

DESCRIPTION 

1-13 

STITLe 

A1Q 

This card is used to label the program output and identify the problem. 





Any number of TITLE cards may be input. Columns 1-30 are printed. 


Card 1. 2-Program Checkout Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SCHECk 

A10 

This card is not recommended for normal production runs. When used, 
the program tries to continue even after fatal errors are encountered. 


Card Set 2.0-Sort Instruction Card 

Note: Card set 2.0 must be input at least once before card set 7.0, $ANAL, the 

analysis card. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SSORT 

A10 

This card directs L225 (TEV126) to call SORT, the response function 
sorting overlay. Cards 3.0 through 6.0 will be read, response 
functions will be sorted, and control will be returned to the 
L225 (TEV126) main overlay. 
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Card Set 3.0-Size and Characteristics of Input Data 
Card set 3.0 must follow card set 2.0 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SIZE 

A10 

Keyword for the size of input. 

11-15 

NRMAX 

15 

When the response functions are to be read from cards, NRMAX is the 
number of responses to be input. 

When the response functions are to be read from a magnetic file, 

NRMAX is the larger of the following two numbers: NDOF, the number 

of generalized coordinates (q's) on the file, and NLD, the number of 
loads on the file. j 

(0 < NRMAX < 100) 

16-20 

NFREQ 

IS 

Number of frequencies to be read from cards or magnetic file. 

Note: All frequency response functions must be defined for the ! 

same array of frequencies. 

(0 < NFRE3 < 250) 

21-30 

FACTOR 

E10.0 

The frequencies read from ca^ds or tape will be multiplied by FACTOR. 
FACTOR is used to convert the units of the input frequencies to Hertz. 
Default is 1.0. 

(FACTOR > 0) 


Card 3.1-Program Checkout Option (Optional) 
If not used, ICKPT1 will be set to zero. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CHECKout 

A10 

This keyword is used for checkout of the SORT program. 

This card is used to print the input matrix data, or the data 
being passed between overlays. 

11-15 

ICKPT1 

15 

= 0 for normal production run (Default) 

= 1 will print the frequency response functions input to the 

SORT routine. 

= 2 will print the input frequency response functions, output 

frequencies and output frequency response functions from SORT. 
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If the frequencies and frequency response functions are to be input on cards, card set 3.0 
(and card 3.1 if used) will be followed by cards 4.0 through 4.4. If they are being input 
by magnetic file, cards 3.0 and 3.1 will be followed by cards 5.0 through 5.2. 


Card Set 4.0-Frequencies and Frequency Response Functions Input From Cards 

Omit cards 4.0 through 4.4 if frequencies and frequency response functions are on 
magnetic file. 


1 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FREOUencies 

A10 

Keyword introducing the card set defining the input frequency array. 

Card 4.1 -Frequency Array 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FI (1) 

7E1C.0 

Elements of the frequency array. 

11-20 

FI(2) 

• 

• 

• 

FI(NFREQ) 


FI ( i ) , i = 1, NFREQ 

Note: F I ( i ) > Fl(i-l) > 0. 

Repeat the card with 7 numbers per card until all elements are 
defined. 


Card 4.2-Response Functions 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

RESPOnse 

A10 

Keyword introducing the response functions being input on cards. 


Card 4.3-Frequency Response Function Array 
Repeat cards 4.3 and 4.4 in pairs NRMAX times. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FRFR(l) 

7E10.0 

Real part of the frequency response function array. 

FRFR(i), 1 = 1, NFREQ 

Repeat the card with 7 numbers per card until all elements are 
defined. 

11-20 

FRFR(2 ) 

• 

• 

• 

FRFR(NFREQ) 
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Card 4.4-Frequency Response Function Array 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

’ 

DESCRIPTION 

1-10 

FRFI(l) 

7E10.0 

Imaginary part of the frequency response function array 
FRFI (i ) , i = 1, NFREQ 

Repeat the card with 7 numbers per card until all elements are 
defined. 

11-20 

FRFI (2) 

• 

• 

• 

FRFI (NFREQ) 


Card Set 5.0-Frequencies and Frequency Response Functions Input From 
Magnetic File 

Omit cards 5.0 through 5.2 if frequencies and frequency response functions are 
input on cards. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FILE 

A10 

Keyword indicating the frequency array and responses are on a 
magnetic file. 

m 

NTAPE 

A10 

File name on which frequency array is found. 

21-25 

IS0L 

15 

The number of the solution on NTAPE which is to be sorted. 
Default = 1. 

26-30 

NSAVEQ 

15 

The number of generalized coordinates (q's) to be saved for 
analysis . 

Note: NSAVEQ < NRMAX. 

Default: All generalized coordinates will be used in analysis. 

31-35 



NSAVEL 

15 

The number of load frequency response functions (L's) to be saved 
for analysis. 

Default: All load responses will be used in the analysis. 

Note: NSAVEL < NRMAX. 
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Card 5.1-Choose q's To Be Saved 

Omit card 5.1 if NSAVEQ = 0 on card 5.0. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

ISAVEQ(l) 

1415 

An array of numbers indicating which q's are to be saved for analysis. 
ISAVEQ(i), i = 1 , NSAVEQ 

Repeat the card with 14 numbers per card until all elements are 
defi ned . 

6-10 

ISAVEQ (2) 

• 

• 

• 

ISAVEQ(NSAVEQ) 


Card 5.2-Choose L’s To Be Saved 

Omit card 5.2 if NSAVEL = 0 on card 5.0. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

ISAVEL(l) 

1415 

An array of numbers indicating which loads are to be saved f cr 
analysis. 

I SAVEL ( i ) , i = 1, NSAVEL 

Repeat the card with 14 numbers per card until all elements are 
defined. 

6-10 



ISAVEL (2) 

• 

• 

• 

I SAVEL ( ilSAVEL ) 


Card Set 6.0-End of SORT Overlay 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1 - 1C 

$END 

A10 

Keyword to indicate the end of data for the SORT overlay. Execution 
is returned to L225 (TEV126) main overlay. 


Card Set 7.0-Analysis Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

SANALysis 

A10 

This card directs L225 (TEV126) to call ANAL, the time history 
calculation overlay. Card *ets 8.0 through 1 1.0 will be read, time 
histories calculated and printed, and execution will be returned 
to L225 (TEV126) main overlay. 
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Card set 8.0-Time Functions 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FUNCTion 

A10 

Keyword indicating the type of time function to be used in the 




analysis. 


11-15 

TYPEX 

A5 

Keyword indicating the type of time function to use. First 5 




characters of Keyword are checked. This includes blank spaces. 




Where: 

STEP function 





TRIPLe step 





IMPULse 




TYPEX = < 

SINE wave 
ONE Minus cosine 
THREE one minus cos. 
RAMPS 





NON-Standard 




See Figure 18 for an illustration of the above time functions. 


Note: Card 8.0 must be followed by cards 8.1 through 8.3 if 

TYPEX ^NON-Standard or by card 8.4 or 8.5 if TYPEX = NON-Standard. 


Card 8.1 -Time Array 

Omit cards 8.1 through 8.3 if TYPEX = NON-Standard. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TCD 

7E10.0 

Array of times used to define the time function. 

11-20 

T(2) 

• 


T(i), i = 1, 7 

61-70 

T( ” 
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Card 8.2- Amplitude Array 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

AHP(l) 

7E10.0 

Array of amplitudes used to define the time function. 

11-20 

61-70 

AMP (2) 

• 

• 

• 

AMP ( 7 ) 


AMP(i), i = 1, 7 


Card 8.3-Period of Time Functions 


Omit card 8.3 if the time function requested is not SINE wave, ONE Minus, or 
THREE one. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PERIOD 

E10.0 

Period is used to define the periodic time functions. Required 
only if the Time Function name is: SINE wave, ONE Minus, or 

THREE one . 


Card 8.4-Tape Containing the Time Functions 

Omit cards 8.4 through 8.5.2 if TYPEX ^ Non-Standard. 
Omit card 8.4 if time function is on cards. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-iO 

TAPE 

A10 

Keyword indicating the Time Functions are on tape. 

11-20 

NTFTP 

A10 

Name of the tape containing the Time Function. 

21-25 

NFILTF 

15 

The numbers of files to skip on NTFTP before trying to read the 
Time Function. 

26-30 

NMATTF 

15 

The numbers of matrices to skip on NTFTP (after file spacing) 
before trying to read the Time Function. 
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Card 8.5-Cards Containing the Time Function 
Omit card 8.5 if time function is on tape. 

Card 8.5 must be followed by cards 8.5.1 and 8.5.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CARDS 

A10 

Keyword indicating the Time Functions are on cards. 

11-15 

NT I ME S 

15 

Number of points in the Time and Amplitude arrays 




( < 1000) 


Card 8.5.1-Time Function Array 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TFT ( 1 ) 

7E10.0 

An array of times used to define the Time Function 

11-20 

TFT ( 2 ) 

• 


TFT(i), i = 1, NTIMES 


• 

TFT (NTIMES) 


Repeat the card with 7 numbers per card until all elements are 





defined. 


Card 8.5.2-Amplitude Array 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TFA(l) 

7E10.C 

An array of amplitudes used to define the Time Function 

11-20 

TFA ( 2 ) 

• 


TFA(i), i = 1, NTIMES 


• 

TFA (NTIMES) 


Repeat the card with 7 numbers per card until all elements are 




defined. 



Card Set 9. 0-Ana lysis Options 


! KEYWORD/ 
! COLS. ! VARIAB1 F. 

FORMAT 

j 

: DESCRIPTION 

, A" iv 

OPT I On s 

.45. SX 

1 Keyword designating the card set containing the analysis options. 

i 11- IS 

I FT ION 

__ 

15 

' See complete description and possible values for IFTION shown in 




; table 3. 

; 16-20 

M 

15 

M = 6, 7, 8, 9 or 10 is used to define the number of points (N) 




1 'n suoroutine HARM (N = 2^), only if 1PTI0N / 1. 

! 



j See next page for more description. 

: ■ 
21-30 

TMAX 

f.10.0 

j Maximum time for which time response is calculated 

! 



' (only for IPTIOli = 3, 6, 7, 8) 

| 31-40 

DELTA! 

E10.Q 

' Time increment (only if I PT I ON * 4, 5, 7, 8) 

1 

j 41-50 

1 

s 

El 0.0 

\ Maximum frequency for which the frequency response functions are 

| 

; 



j considered to be defined. 


A summary of the available options is snown in table 3. A discussion of the options is 
presented in section 5-3. 


Card 9.1-Conjugate Response Function (Optional) 



! KEYWORD/ | I 

i ..Ols. 1 Vm£ 1A31E FORMAT j 

i i I 

j 1-10 j CONJUg atc j A10 



_ DESCRIPTION 

Keyword to cause the input frequency response functions to be 
conjugated before transform is performed. 

Default is no conjugation. 




Table 3.— Summary of Ave 


Option 

Input data 

1 

None 

2 

M 

3 

M, t 
1 MAX 

4 

M,At 

5 

M.At 

6 

M ’ T MAX 

7 

M. T MAX' At 

8 

M ' T MAX' At 


Parameter Options 


Restrictive conditions 

F HARM =F MAX1 = F MAX0 
M(1) = 10 

T MAX =1024/F MAX0 
At=1/F MAX0 
F HARM =F MAX1 = F MAX0 

t max =n/f maxo 
A t = 1/f maxo 

f harm =f maxi <f maxo 
t max > n/f maxo 

At>1/F MAX0 

F HARM =F MAX1 <F MAX0 

t max >n 'f max0 

At>1/p 

h MAX0 

F HARM =F MAX1 =F MAX0 
IF T MAX <N/F MAX0 (usual si,uation) 
At<N/NN*F MAXQ 
or 

F HARM =F MAX1 <F MAX0 

,F T MAX >N/F MAX0 
At>N/NN*F MAX0 

F HARM >F MAX1 

f maxo >f maxi 

At is changed only as a last 
resort 

F HARM >F MAX1 

F MAX0 >F MAX1 

Tmax ' s changed only as a last 

resort 
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Card 9.2-Integration Selector (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

_ _ - - 

DESCRIPTION 

1-10 

INTEGrate 

A10 

Keyword 

introducing the integration selector. 

11-20 

TYPEX 

A10 

The type of integration to be performed in takinq the inverse j 




i 

transform where TYPEX must be one of the following keywords. 1 





COMPLex (Default) 




TYPEX =. 

REAL 





IMAGInary 


KEY WORD 

EQUATION USEO 

COMPLex 

r 2rFMAXl_ 

x(t) = 2 Real ~ f x(u>)e 1ut du> 

^ 0 ' 

IMAGInary 

2tt FMAX1 

x (t) = ~hr ^ ma 9 1 | sir 

0 

REAL 

2ttFMAX1 

x (t) = 4 j *Real | x(w) ^ cos cot dw 

0 


Card 9.3-Print Parameter Changes 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 1 DESCRIPTION 

! 

1 1-10 

PRINT 

— 

A10 

. 

Keyword requesting printing of control data changes made by the 
program to TMAX or DELTAT. 

Default is not printing TMAX or DELTAT changes. 
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Card Set 10.0-Output Specifications 


jCOLS . 

1-10 

11-20 

21-25 


KEYWORD/ 

VARIABLE 


FORMAT 


DESCRIPTION 


OUTPUt 


A10 


STANDa rd or A1 0 
NON-St andard 

ICKPT2 15 


Keyword to specify output options. 

If keyword STANDARD is used the only output printed is 
Time Histories. Omit Cards 10.1 through 10.4. 

If keyword NON-STANDARD is used the desired options must be 
specified on Cards 10.1 through 10.4. 

This variable is used for check out purposes in the program. 

= O for normal production run. 

= 1 will print input time functions from the analysis routine. 

= 2 will print input frequencies and input responses from the 

analysis routine. 


Card 10.1-Print Options 

Omit cards 10.1 through 10.4 if STANDard keyword is used on card 10.0. 

Card 10.1 may be repeated with different keywords if more than one nonstandard 
output type is desired. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PRINT 

A10 

Keyword for print options. 

11-20 

! 

TYPEX 

A10 

Keyword to indicate the desired print option. TYPEX must be one 

of the following: 

TIME - time histories of the desired output functions with the time 
history of the forcing function printed with the first output 
quantity only. 

INTERpolated - print interpolated transfer functions used in 
convolution. 

ALL Transfer functions - print interpolated transfer functions, 

transformed forcing function and product of interpolated functions 
and forcing function. 

Note: ALL Transfer functions option includes INTERpolated option. 

MAXIMum - tninimums - print maximum and minimum values of output 
functions along with times at which they occur. 
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Card 10.2-Selected Output Times (Optional) 


Omit cards 10.2 and 10.2.1 if the time response is to be output at the times at 
which it was determined in the analysis. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SELECt 

A10 

Keyword indicating time responses to be printed at selected 
output times as specified by user. 

11-15 

NTPTS 

15 

Numbers of time points at which the user chooses time responses 
to be printed 
(NTPTS < 1024) 

16-25 

DELTO 



E10.0 

Value of At which user specifies along whith NTPTS to determine 
the output time array. 

TJ(i) = i * DELTO 
where i ■ 1, NTPTS 


If DELTO = 0.0, user must read in time array on card 10.2.1. 
Card 10.2.1-Time Array for Output 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TJ(1) 

7E10.0 

Time array for output of time response. 

11-20 

TJ(2) 

• 

• 


TJ(i), i = 1, NTPTS 


• 

TJ(NTPTS) 


Repeat the card with 7 numbers per card until all elements are 




defined. 


Card 10.3-Plot Time History (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PLOT 

A10 

Keywords to indicate that the time histories are to be saved for 
possible plotting. (See tape PLTRES in Section 6.6.2). 

11-20 

TIME 

A10 
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Card 10.4-Plot Function (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-15 

PLOT 

A10 

Keywords to indicate that the output frequency response functions are 
to be saved for possible plotting. (See tape PLTRAN in Section 6.6.2). 

11-20 

TRANS 

A10 



Card Set 11.0-End of Analysis 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

A10 

Keyword to indicate the end of data for the ANAL overlay. 
Execution is returned to L225 (TEV126) main overlay. 

Card Set 12.0-End of Program 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SOU IT 

A10 

This is the last card input to L225 (TEV126). 




It terminates the execution of the program. 
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Requirements 
or Function 


<~n 

4s. 





$T IME 
tTITLe 
SCHECk 
$S0RT 



Input Data Size and 
Characteristics 

SIZE 

Program Checkout 
(Optional ) 

CHECKout 




Freq. Input on Cards 

FREQUenci es 

Frequency Array 

FI ( 1 ) 

Response Functions 
Input on Cards 

RESPOnse 

1 

Real Part of 
Response Functions 

FRFR(l) 

Imaginary Part of 
Response Functions 

FRFI(l) 



Magnetic File Options 

FILE 

NSAVEOO; CS 5.0 

ISAVEQO ) 


Keywords and/or Variables 



6.5.3 SUMMARY OF CARD INPUT DATA 













Requirements 
or Function 

T" 

Keywords and/or Variables 

Card 

Format 

Reference 
Card Set (CS) 


Begin Input of Analysis Parameters 




SANALysis 

A10 

7.0 

Selection of 
Time Function 

FUNCTion TYPEX 

A10.A5 

8.0 


Omit Cards 8.1 through 8.3 if 
TYPEX = NON-Standard; CS 8.0 



Time Array 

T(l) T(2) • • • T( 7) 

7E10.0 

8.1 

Amplitude Array 

AMP { 1 ) AMP ( 2 ) • • • AMP{7) 

7E10.0 

8.2 


Omit Card 8.3 if TYPEX / SINE wave, ONE Minus, 
or THREE one; CS 8.0 



Feriod of Time 
Function 

PERIOD 

. 

E10.0 

8.3 


Omit Cards 8.4 throuqh 8.5.2 if 
TYPEX f NON-Standard; CS 8.0 




Omit Card 8.4 if time function is 
input on cards 



Time Function on 
Magnetic File 

TAPE NTFTP NFILTF NMATTF 

2A10.3I5 

8.4 ' 


Omit Cards 8.5 through 8.5.2 if 
time function is input on magnetic file 



Time Function 
on Cards 

CARDS NT I ME S 

A10.I5 

8.5 

Array of Times 

TFT(l) TFT(2) • • ■ TFT(NTIMES) 

7E10.0 

8.5.1 

Array of Amplitudes 

TFA(l) TFA{2) • • • TFA(NTIMES) 

7E10.0 

8.5.2 




Requirements 
or Function 

Keywords and/or Variables 

Card 

Format 

Reference 
Card Set (CS) 

Analysis Options 

OPTIOns I PT ION M TMAX DELTAT FMAXO 

A10.2I5 

3E10.0 

9.0 

Conjugate Response 
Functions (Optional) 

CON JUgate 

A10 

9.1 

Integration Selector 
(Optional ) 

INTEGrate TYPEX 

2A10 

9.2 

Print Parameter 
Changes (Optional) 

P.RINT 

A10 

9.3 

Output Specifications 

OUTPUt STANDard or ICKPT2 

NON-Standard 

2A10.I5 

10.0 


Omit Cards 10.1 throuqh 10.4 if keyword STANDard 
appears on Card 10.0. Card 10.1 may be repeated 
with different keywords if more than one non-standard 
output type is desired 



Print Options 

PRINT TYPEX 

2A10 

10.1 


Omit Cards 10.2 and 10.2.1 if time response is to be 
output at the times at which it was determined in 
the analysis 



Output Time Selection 
(Optional ) 

SELECt NTPTS DELTO 

A10 ,15 
E10.0 

10.2 

Output Time Array 
DELT0=0; CS 10.2 

TJ(1 ) TJ(2) • • • TJ(NTPTS) 

7E10.0 

10.2.1 

Plot Time Histories 
(Optional ) 

P.LOT TIME 

2A10 

10.3 

Plot Output Response 
Functions (Optional) 

PLOT TRANS 

2A10 

10.4 

End of SANALysis 
Overlay 

SEND 

A10 

11.0 

End of Program 

SQU I T 

A10 

12.0 



6.5.4 MAGNETIC FILES (TAPE OR DISK) 


Optionally, program L225 (TEV126) can read the frequencies and responses it processes 
from a magnetic file "NTAPE”. Card set 5.0 (see sec. 6.5.2) is used to specify whether 
the responses are on cards or on "NTAPE”, which is described in figure 19. 

L225 (TEV126) may also read time functions from the magnetic file "NTFTP”. Card 8.4 
is used to indicate that the time functions are on "NTFTP”. See figure 20 for a 
description of' NTFTP”. 

Both of the files must be written in the READTP/WRTETP format. 1 


"NTAPE”-Frequency Response Input Tape 

"NTAPE” is a magnetic file (tape or disk) containing the frequency-dependent responses 
to be analyzed by program L225 (TEV126). "NTAPE” must be written in the 
READTP/WRTETP format. The file is usually generated by the program L221 (TEV156) 
of reference 5. "NTAPE” may contain several logical files, one for each solution on a 
"STANDARD” tape generated by L221 (TEV156). The contents of a typical logical file 
are displayed in figure 19. 


"NTFTP”-Time Function Input Tape 

"NTFTP” is a magnetic file (disk or tape) containing the time function used in the 
convolution with the frequency response functions. "NTFTP” is specified via card 8.4. 
The file must be written in the READTP/WRTETP format. 

"NTFTP” may contain several logical files and matrices. The file and matrix skipping 
parameters on card 8.4 enable the program to read the two required matrices, the array 
of times, and the array of amplitudes. See figure 20 for the contents of "NTFTP”. 


1 Clemmons, R. E.: Programming Specifications for Modules of the Dynamic-Loads System 
to Interface With FLEXSTAB. NASA contract NAS1-13918; BCS-G0701, September 1975. 
(Internal document) 



SIZE 


Where : 


1 x NFREQ 
1 x NDOF 
1 x NDOF 
1 x NLD 
1 x NLD 
1 x NDOF 
1 x NDOF 
1 x NLD 
1 x NLD 


1 x NDOF 
1 x NDOF 
1 x NLD 
1 x NLD 



NFREQ = no . of frequencies 

NDOF = no. of degrees 
freedom 

NLD = no. of loads 

{FREQ} = frequency array in 
cycles per second 

{Qd} and .{Qi} are the real 
and imaginary frequency 
response function arrays. 

|Lp} and {Lj} are the real 
ana imaginary load frequency 
response function arrays. 


Figure 19. — "NT APE"— Frequency and Response Input Tape 
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SIZE 
1 X IN 
1 x IN 


Preceded by 

NFILTF logical files 

and NMATTF matrices 



Where : 

{ TFT } = gust time array 

{ TFA} = gust time function 
array 

IN S 1024 

The elements of {tFt} must 
be monotonically increasing 
and TFAi must correspond 
with TFTi 


NFILTF and NMATTF are 
described in Section 3.5.1 
under Card 8.4 


Figure 20.-"NTFTP"-Time Function Input Tape 
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6.6 OUTPUT DATA 


L225 (TEV126) generates printed output and, optionally, magnetic files (disk or tape). 
6.6.1 PRINTED OUTPUT 


Each time L225 (TEV126) is executed, the following items will be printed: 

• A one page block indicating the program name and version plus the current date 
and time 

• All card input data plus interpretive comments 

• Execution control data generated by the program 

• Solution frequencies after being scaled according to card input 

• Time histories 

• A one page block indicating that the program execution has ended. 

In addition, the following items may be printed if requested by card input. 

• Changes made to the program control parameters TMAX and DELTAT (see 
card 9.3) 

• Interpolated frequency response functions used in the convolutions (see card 10.1) 

• All frequency functions: the interpolated frequency response functions, transformed 
input time functions, interpolated output frequency response functions, and 
interpolated input time functions (see card 10.1) 

• The times at which minimums and maximums are attained 

For checkout purposes, two other print options are provided: 

• Card 3.1 requests the printing of intermediate data in SORT-the sorting overlay 
(L225,l,0) 

• ICKPT2 of card set 10.0 requests the printing of intermediate calculations in 
ANAL-the analysis overlay (L225,2,0) 
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6.6.2 MAGNETIC FILES (TAPE OR DISK) 


Optionally, L225 (TEV126) will write two magnetic files containing plot data; 
"PLTRES” and "PLTRAN”. 


Card 10.3 is used to request the preparation of the file "PLTRES” containing gust time 
functions and output time responses (see fig. 21). 


Card 10.4 is used to request the preparation of the file "PLTRAN” containing 
transforms of the output convolutions (see fig. 22). The two files are both written in the 
READTP/WRTETP format. 


Both the gust time functions and the output time responses are written on file 
"PLTRES” if requested by card 10.3. "PLTRES” will contain alternate arrays of 
independent and dependent variables suitable for plotting. 


The absolute values of the transform of the output convolutions are written on 
"PLTRAN” if requested by card 10.4. "PLTRAN” will contain alternate arrays of 
independent and dependent variables suitable for plotting. 


6.7 RESTRICTIONS 


The limitations placed on input variables are described in detail with each data card in 

section 6.5.2. The primary restrictions are summarized in the following listing. 

• SORT overlay must be executed before execution of the overlay ANAL. 

• Maximum number of input response functions is 100. 

• Maximum number of times and amplitudes used to define a standard time function 
is 100. 

• "NON-STANDARD” type time function may be defined with up to 1000 time 
points. 

• Maximum number of time points at which the time histories can be printed 
is 1000. 
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SIZE 


Where 


1 x MM 

{Tlj 

1 x MM 

M 

1 x MM 


TI 

or 


or 

1 x NTPTS 

1 

TJ 

1 x MM 


or 





l Convolution 
[ No. 1 


1 x NTPTS 


TI = time array 

FT = time function array 

TJ = time array 

X = convolution time 
response 

MM £ 1024 
NTPTS £ 1024 


1 x MM 


TI 


or 


or 


All output files (other 
than print) ire written 
with subroutine WRTETP. 


1 x NTPTS 
1 x MM 
or 

1 x NTPTS 


TJ 


> Repeated for 
I each Convolution 


NOTE: This same file structure is 

repeated for each analysis if 
PLOT TIME is requested on 
Card 10.3. 



Figure 21. ~"PLTRES"-File of Time Responses for Plotting 
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SIZE 
1 x MM 
1 x MM 
1 x MM 
1 x MM 


1 x MM 
1 x MM 


Where : 



Convolution 
No. 1 


Convo lut ion 
No. 2 


F = frequency array 

ABSX = transform of the 
convolution 


MM < 1024 


All output files (other 
than printed) are written 
with the subroutine WRTETP . 


Repeated for 
each convolution 


NOTE : This same file structure is repeated for each analysis if 

PLOT TRANS is requested on Card Set 10.4. 


Figure 22.—"PLTRAN”—FHe of Response Frequency Functions for Plotting 
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6.8 DIAGNOSTICS 


All errors detected by the program L225 (TEV126) will cause the program to print a 
diagnostic message. These messages are self-explanatory. The first line printed is in the 
format: 

******** FATAL ERROR nn name 

where nn is the error number (from 1 to 17), and name is the name of the routine in 
which the error was discovered. 

The following list is a brief description of each error: 

1 Premature end-of-file was encountered in the card input file. 

2 Keyword is not appropriate. 

3 NFREQ (number of frequencies) must be between 1 and 250. 

4 NRMAX (maximum number of responses) must be between 1 and 100. 

5 FACTOR (frequency scaling factor) must be greater than zero. 

6 Data or option falls outside the permissible range. 

7 ISAVEQ(i) (array of Q’s to be saved) and ISAVEL(i) (array of loads to be saved) 
must fall inside the range: 

0<ISAVEQ(i) NRMAX 
0<ISAVEL(i) ^ NRMAX 

8 FI(i) (array of frequencies) must be increasing in magnitude FI(i+l) > FI(i) 

9 Error returned from subroutine FETAD, error with file "NTAPE”. 

10 Error returned from subroutine READTP (see READTP error codes which follow). 

11 Matrix read is not of the expected size. 

12 Error returned from subroutine WRTETP (see WRTETP error codes which follow). 

13 Field length too short. 

14 Error returned from subroutine INITAL * . 

15 Error returned from subroutine TFRESP *. 

16 Error returned from subroutine CONVOL * 

17 If FHARM (maximum frequency) is nonzero, then DELTAT (time increment) must 
be zero, or vice versa. 
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6.8.1 READTP ERROR CODES 


'I 


IRR - 0 No errors were detected. 

= 1000+i End-of-information was encountered after skipping (NFILES-i) files 
of the "NFILES” files to be skipped. 

= 2 Number of matrices to be skipped (NMATS) was less than zero. 

= 3000+i End-of-file was encountered after skipping (2 *NMATS-i) records of 

the (2 *NMATS) records to be skipped. 

= 4 NROWS > NROWD, where NROWS is the actual number of rows in 

the matrix being read, and NROWD is the dimensioned row Bize of 
the array in core. 

= 5 Name check failed. 

= 6 Either NROWS * NCOLS = 0 or NROWS *NCOLS is greater than 

the allowable buffer size. NCOLS is the column size of the matrix 
being read. 

= 7 End-of-file was read where the matrix was expected. 

6.8.2 WRTETP ERROR CODES 

IRR =0 No error was detected. 

= 1000+i End-of- information was encountered after skipping (NFILES-i) files 
of the "NFILES” files to be skipped. 

= 2 Number of matrices to be skipped (NMATS) was less than zero. 

= 3000+i End-of file was encountered after skipping (2 *NMATS-i) records of 

the (2 *NMATS) to be skipped. 

= 4 Either NROWS *NCOLS = 0 or NROWS * NCOLS is greater than 

the allowable buffer size. 
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7.0 SAMPLE PROBLEM 


7.1 PROBLEM DESCRIPTION 

The sample problem presented in this section is a symmetric analysis of a large 
transport airplane. Two rigid body degrees of freedom, vertical translation and pitch, 
plus six elastic degrees of freedom constitute the eight generalized coordinates. Eleven 
load responses were considered, three wing shears, three wing bending moments, four 
wing accelerations, and the airplane acceleration at the center of gravity. Retained for 
the sample problem were the pitch and first elastic modal response (Q’s 2 and 3), and a 
wing acceleration response (L 10). Frequency response functions were calculated at 210 
frequency values. The airplane was subjected to a one-minus cosine wave of unit 
amplitude, having a period of 0.82 sec. 
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^sm 


(TIME HISTURY ANALYSIS 
(TITLE 

(TITLE L225ITEV126I OUCUHENT SAMPLE PROBLEM. 

(TITLE 

(SORT FREQUENCY RESPONSE FUNCTIONS 

(TITLE STANDARD INPUT TAPE FROM L22UTEV156) CONTAINS 19 FREQUENCY 

(TITLE RESPONSE FUNCTIONS (8 GEN. COORD. AND 11 LOADSI . 

(TITLE THE 2ND AND 3RD COORDINATE RESPONSES 

(TITLE AND THE 10TH LOAD RESPONSE MILL BE RETAINED. 

SUE 11 21Q 1.0 

FILE IRTAPE 1 2 l 

2 3 

10 

(END SORT OF FREQUENCY RESPONSE FUNCTIONS 
(ANALYSIS OF TRANSFER FUNCTIONS 
FUNCTION SINE WAVE 
0.0 .82 
1.0 
.82 

OPTIONS 6 103.0 0.0 0.0 15.0 

INTEGRATE COMPLEX 

PRINT CHANGES IN CONTROL DATA 

OUTPUT NON-STAND. 

PRINT TIME HISTORIES 
PRINT ALL TRANSFER FUNCTIONS 
PRINT MAXIHUMS AND MINIHUMS 

SELECT 100.03 

(END ANALYSIS 

(QUIT TIME HISTUKY PROGRAM 


ON 

■-J 


7.2 LISTING OF SAMPLE PROBLEM INPUT DECK 



*»4 

******** 

prugram 

L 225b 3 

VERSION 76/09/20 



BEGINNING EXECUTION ON THE BCS O600. 



CATE OF 

RUN IS 

76/11/03. 


* + 1 

TIME OF 

RUN IS 

»»**•*** 

12.32.26. 



I IT [ ME HISTORY ANALYSIS 
( (TITLE 

(ST1TLE L225ITEV126I DOCUMENT SAMPLE PROBLEM. 
( (TITLE 

(SSOKT FREQUENCY RESPONSE FUNCTIONS 



7.3 LISTING OF SAMPLE PROBLEM PRINTED OUTPUT 



SORT HAS BEEN ENTERED TO SORT FREQUENCY 
RESPONSES TO BE ANALYSED BY ANAL • 


UTITLE STANDARD INPUT TAPE FROM L221 ( TE V l 5b) CONTAINS 19 FREQUENCY 

ItTITLE RESPONSE FUNCTIONS 18 GEN. COCRD. AND 11 LOADS! . 

CiTITLE THE 2ND AND 3RD COORDINATE RESPONSES 

((TITLE AND THE IOTH LOAD RESPONSE HILL BE RETAINED. 

(SIZE 11 210 1.0 

THE MAXIMUM NUMBER OF INPUT RESPONSES IS II 
THE NUMBER OF INPUT FREQUENCIES IS 210 

THE FREQUENCY SCALING FACTOR IS 1.0000 


FILE IRTAPE 121 

THE RESPONSES HILL BE READ FROM IRTAPE 


THE 

SOLUT ION 

TO BE 

SORT 60 

IS 

NUMBER 1 

THE 

2 Q-S TO 

BE RETAINED 

FOR 

ANALYSIS ARE 


2 3 





THE 

1 LOADS 

TO BE 

RETAINED FOR ANALYSIS ARE 


10 





THE 

CURRENT 

FIELO 

LENGTH 

IS 

115000 

THE 

REQUIRED 

FIELD 

LENGTH 

IS 

055421 



o 


THE FREQUENCIES WRITTEN UN SRTFRP AFTER BEING MULTIPLIED BY A FACTOR OF 
FREQUENCIES ARE IN UNITS OF HER T l 


7958E-02 

.1592E-GI 

.3183E-Q1 

.*7756-01 

•0 36oE-Gl 

•7162E-01 

.7958E-01 

.875*6-01 

•95*9E-01 

. 1035E + 00 

aiuE^oo 

.1273E+0Q 

.1*326+00 

• l 59ZE+00 

. L751E+0Q 

19 1 OE 4-00 

.20696*00 

.22286+00 

• 2 387E+00 

• 25*66+Gu 

2706E+00 

.28656*00 

• 302*E+Q0 

• 3 1 83 E +00 

•33*26+00 

•3501E+00 

.38206+00 

•*297E+00 

.*7756+00 

.52526+00 

•57306+00 

• 6 2076+00 

.66856+00 

.71626+00 

•76396*00 

•79586 +00 

.83566+00 

•875*E+00 

.9151E+00 

•95*9E+0G 

99*76+00 

• 1 035E+0 1 

.107*6+01 

.lll*t+Ol 

.115*6+01 

1 19*6+01 

• 1 233E +01 

• 12 73E+01 

.13136+01 

• 1353E+01 

1393E+01 

• 1 *32E +01 

•1*726+01 

• 1 5126+01 

.15526+01 

.15926+01 

• 1631E + 0L 

. 16716+01 

.17116+01 

• 1751E+01 

1790E+01 

.103OE+OI 

.18706+01 

.19106+01 

. 19506+01 

19896+01 

•2029E+01 

.2069E+01 

.21096+01 

.21*96+01 

■ 2 188E+QL 

.22286+01 

. 22 68E + 01 

•2308E+01 

•23*86+01 

.23876+01 

• 2*276+01 

• 2*67E+QL 

•2507E+QL 

.25*66+01 

.25866+01 

.26266+01 

• 2666E+Q1 

.2 7066+01 

.27*56+01 

•2785E+01 

.2025E+O1 

• 28656 + 01 

•2905E+01 

.29**6+01 

• 298*E+0 1 

.302*6+01 

• 306*E+0 l 

•310*6+01 

•31*36+01 

•3183E+01 

.3223E+01 

.32636+01 

.33026+01 

.33*26+01 

•3382E+01 

.3*226+01 

•3*62E+01 

•3501E+01 

.35*16+01 

.358 iE+O l 

•36216+01 

.36616+01 

•3700E+01 

.37*06+01 

3780E+01 

.38206+01 

.36606+01 

•3899E+01 

.39 396+01 

• 3979E + 0 1 

• *0t 9E +01 

.*0586+01 

•4098E+01 

.*1386+01 

• *1 78E+01 

.*2186+01 

. *257E+0l 

•*297t+01 

.*3376+01 

• *377E *0 l 

.**176+01 

• **56E + 01 

•**96E+01 

.*5366+01 

.*5766+01 

.*6156+01 

• *655E+0l 

•*695E+01 

.*7356+01 

•*7756+01 

•*ai*E+oi 

•*85*£+0l 

• *89*E +0 1 

.*93*6+01 

.*97*6+01 

.50136+01 

.50536+01 

.50936+01 

.51336+01 

• 5173E+01 

•5212E+01 

•5252E+01 

.52926+01 

.53326+01 

.537 16 + 0 1 

« 5*1 IE +0 1 

.5*516+01 

• 5*9 IE +01 

.55316+01 

•55706+01 

• 56506+01 

.57306+01 

.58096+01 

•58896+01 

.59686+01 

•60*06+01 

• 61 27E + 01 

•6207E+01 

•6267E+01 

.63666+01 

•6**6E +01 

.65256+01 

. 6605E +0 1 

•66056+01 

.676*6+01 

• 68**E +01 

.69236+01 

.70036+01 

.70826+01 

.71626+01 

.72*26+01 

• 7321E+01 

.7*016+01 

.7*806+01 

.7560E+01 

.76396+01 

• 77 19E+01 

• 7799E+01 

• 78 78 E+Ol 

.79586+01 

• 8037E+01 

.81 17E+01 

.81966+01 

.82766+01 

. 8356E *0 1 

.8*356+01 

.85156+01 

.859*6+01 

.067*6+01 

.875*6+01 

.88336+01 

. 8913E+01 

.89926+01 

• 90 72E + 01 

•9151E+01 

•9231E+0L 

•931 1E+01 

.99*76+01 

.10986+02 

. L194E+02 

.13056+02 

. l*01E+02 

.1*966+02 

.15926+02 


(JENU SORT CF FREQUENCY RESPONSE FUNCTIONS 
1 i ANAL YS IS OF TRANSFER FUNCTIONS 


1.0000 



ANAL HAS BEEN ENTERED TO ANALYSE 
FREQUENCIES AND RESPONSES FROM SORT 


(FUNCTION SINE WAVE 
THE TIME ARRAY IS 

0. • 8200E+00 0. 0. 0. 

THE AMPLITUDE ARRAY IS 

• 1000E >01 0. 0. 0. 0. 

PERIOD OF TIME FUNCTIONS IS .82Q0E+00 

I OPT I ONS 6 103.0 0.0 0.0 15.0 

THE SELECTED ANALYSIS OPTIONS ARE 
OPTION SWITCH (IPTIONT - 6 

(HI « 10 

MAXIMUM T (ME RESPONSE (TNAXI * .3000E*0l 

TIME INCREMENT (OELTATI * 0. 

MAXIMUM FREQUENCY I FHARM ) - 0. 

MAXIMUM FREQUENCY (FMAXO) - .1500E+02 

( INTEGRATE COMPLEX 

(PRINT CHANGES IN CONTROL OATA 

(OUTPUT NON-STANO. 

(PRINT TIME HISTORIES 

I PRINT ALL TRANSFER FUNCTIONS 

(PRINT MAX (MUMS ANO HlHIMUHS 

I SELECT 100.03 

( SENO ANALYSIS 

SUBROUTINE EQUAL 

TMAX CHANGED FROM .300000E*01 TO 




0 . 0 . 

0 . 0 . 

i 


299665E+0 & . 



-0 

to 


EXECUTION CUNIKUL OATA GENERATEO BY THE PROGRAM 


FREQUENCY INCREMENT IUELTAFJ * 
MAX. FREQUENCY IFMAXt) « 
MAX. FFT FREQ • ( FHARH ) = 
MAX. INPUT FREQ. ( FMAXO) » 
TIME INCREMENT (DELTATJ « 


NO. OF INPUT FREQ (NFREQ) » 210 
NO. OF TIME INTERVALS 4 NN ) » 221 
NO. OF TIME POINTS 4 MM 1 * 222 
ELEMENTS IN ARRAY X IN) « 1024 
ELEMENTS IN ARRAY TFT ( INT « 0 
NO. OF FREQ. POINTS (NPQINT1 « 210 


7201 576E-01 
1591 54BE+02 
73744 14E+02 
1591 548E+02 
1 356 040E-0 1 





% CONVOLUTION NUMBER 


1 * 



I 

FREQUENCY 

INTERPOLATED INPUT FREQUENCY 
RESPONSE FVJNCTl’GN 
FRF 

TRANSFORM OF 
IMPULSE TIME 
T FRF 

THE GUST 
FUNCTION 

TRANSFORM OF THE 
TIME CONVOLUTION 
FRF • TFRF 



REAL IMAGINARY 

REAL 

IMAGINARY 

REAL IMAGINARY 


l 

0. 

-.9221386-03 

-.4828746-06 

-.1755386-03 

0. 

.1616716-06 

.8476306-10 

2 

.7201586-01 

-.91761 0E-Q3 

.2834286-03 

.8740436-02 

• 475066E-01 

-.2 14850 E- 04 

-.4111526-04 

3 

• 1440326*00 

-.6162046-03 

.6764596-03 

•346402E-01 

.894B67E-01 

-• 889865E-04 

-.4953746-04 

4 

• 2160476*00 

-.3220176-03 

.1045046-02 

.7505796-01 

• l 2092 1E*00 

-.1505376-03 

.39499oE-04 

5 

i 2 000636*00 

.2034416-03 

• 8 12 1 36E-03 

.1261376*00 

. 13 7759E *00 

— . 76 L264E-04 

• 1 414876-03 

6 

. 3600796*00 

.3934116-03 

• 40**93 LE-0 3 

. 182966 E*QO 

.1372046*00 

.1639156-04 

. 1281156-03 

7 

.432C95C*0C 

.3212036-03 

• l 04 1 3 7E -03 

.240127C*00 

.1163026*00 

.5533456-04 

• 8 22 3 *»6E— J 4 

8 

• 504 l 10E *00 

.2404286-03 

* B56034E— 04 

.2919926*00 

.6153786-01 

.632007E-04 

.44657oE-04 

9 

• 5 76 1 266*00 

.1 768366-03 

• 3 475 86F -04 

• 333429 C *00 

.2897406-01 

.5795516-04 

. 16 71 32 1-04 

10 

.6401 426 *00 

.13071 LC-03 

.9826206-05' 

.3601606*00 

-.35746GE-01 

.4742626-04 

-.1133776-05 

11 

.7201586*00 

.9464936-04 

- • 44262 8E-05 

.3691536*00 

L07966E+00 

.34462 16-04 

1 18529 fc-04 

12 

. 792 1736*00 

.6577076-04 

-. 1 184246-04 

.3588766*00 

182353E*00 

.2144416-04 

-.1 62434c- 04 

13 

• 0641896*00 

.4041466-04 

r. 1460006-04 

• 3 294 L5 6*00 

- . 2 53-.00E *00 

.9562536-05 

—.151 1 o4E-04 

14 

• 9 162G5E *0 C 

.1545166-04 

-.1310716-04 

.20243 16*00 

-.31 5899t *00 

• 2234 79c- C6 

- • 3 5o255 c- j 5 

15 

. 1000226*01 

-.1144266-04 

-.239^586-05 

.2209926*00 

-.3653846*00 

-.3404396-05 

. 365 L 31 1- J 5 

16 

. 10e024E*Ol 

-.3506936-04 

.297709004 

.1492706*00 

-.3985006*00 

.6628936-05 

• l 64 1 9C 6-u4 

17 

. 1152256*01 

-.1608066-04 

.6470716-04 

.7215456-01 

— .4 1 32 58E *00 

.3384506-04 

• 1 2 76 CSE-0 4 * 

ie 

. 122h27C+01 

.4093726-04 

.3131956-04 

-.5 183586-07 

- • 4C 9 l 59 6 *00 

.3306026-04 

-.1717186-04 

19 

L29620E *01 

.5848056-04 

• 4934626 -04 

-.7770336-0. 

-.3071 82 fc* 00 

.1455636-04 

- . 2 04 759E- J4 

20 

«i36C30E»0 ! 

.5547156-04 

• 2 79t>3 16-04 

-.1411616*00 

— • 34 963uE *00 

• 1 94c5 CE- 05 

-.2334216-04 

21 

i 1440326*01 

.4862926-04 

• 1 612936-04 

1919036*00 

„- *2,999056 *00 

-.4494olE-05 

-. 1 76 7 6-04 

22 

• 1512336*01 

.4203506-04 

.9325016-05 

-.2276896*00 

-.2421006*00 

-.7313246-05 

— . l2Jvuo2t — j4 

£3 

. 1 584 366*01 

.3638076-04 

.52518L6-C5 

-.2475216*00 

-.1607026*00 

-.8055566-05 

- • 7 d 7 i 58t-05 

24 

• lo5b JcE *01 

.3155376-04 

.2624366-05 

-.25 1709 E*U0 

-.1200386*00 

7627186-05 

-. 4449tOE-05 

25 

• 1 72H 386 *0 1 

.2739266-04 

.8937576-06 

-.2417516*00 

-.642398E-01 

-.65644 76-0*5 

197696E-j5 

26 

.1800396*01 

.2374436-04 

-.2219706-06 

-.2201246*00 

- . 1 o 39 l 3E-0 1 

-.5230356-05 

- • 340j 4Cc-06 

27 

.1372416*01 

.2047396-04 

-.8992116-06 

-.1899886*00 

.21 1 8 19 6-0 1 

-.367076E-05 

•604516b- Jo 

28 

. 1944436*0 l 

.174663E-04 

-.1188996-05 

-.1548546*00 

.4720116-01 

-.2640526-05 

. IC0955C-05 

29 

.2016446*01 

. 14596C6-04 

-.1075276-05 

-. 1 L8245E*00 

.6179446-01 

1 0554 5t -05 

. I0*9u9fc-u5 

30 

.208046E* 01 

• L l 7 72 56-04 

-.4047446-06 

-.8337646-01 

.6560416-01 

- * 954994 Ei-06 

« fc Oc>u o5 fc— u 6 

31 

• 2 1604 7E *01 

.5045656-05 

.1224376-05 

5 29006 E-0 1 

.6040866-01 

-.5 524 94e|-96 

• 49lo 77c-06 

32 

•223249E*0l 

.7G83456-C5 

.4529816-05 

-.28 7 165 E-0 L 

.4848126-01 

-.423G23C-06 

. 2 1 33 J4E- 06 

33 

• 2 304 50 1 ♦ 0 1 

.84861 46-05 

.9359296-05 

-.1187466-01 

.3235576-01 

- .4052666-06 

.1633016-06 

34 

.2376526*01 

.1423126-04 

.10325o6-04 

-.2571326-02 

• 1 4 7*«6 76 -0 1 

108 eo 2 E -06 

. 1033136-vo 

35 

.244854E*01 

. 1 73 74 OE -04 

• 63994 5E -05 

-.2275956-03 

- . 2 1 0540E -02 

.9515246-00 

- . 3 603 5cE-37 

36 

. 2 520 55E *01 

.1698126-04 

.2705 706-05 

-.3633066-02 

-.1624176-01 

-.1783336-07 

-.285o47£-o6 

3? 

.2592576*01 

.1547026-04 

.5102716-06 

-.11 1606 E-0 L 

-.2^.32726-01 

1593476-06 

4 129d6fc-0o 

38 

. 2664586*0 l 

.1388786-04 

-.7248656-06 

-.2097866-01 

-.3 L6786E-0L 

-.3143106-06 

- . 42**74LL-jo 

39 

• 2 73ut>CE *0 1 

.1244926-04 

144016E-05 

-.3124166-01 

-.3223706-01 

-. 43538 t«F- 06 

-.35o551 c-06 

40 

.2803ol6*0 1 

. 1 1 1 345E-04 

-. 1 071636-05 

-.4034276-01 

- . 2 35909 6-0 L 

-.5027336-06 

- • 24^9 1 oc-06 

41 

.2800636*01 

.9894246-05 

-.2C863O6-05 

-.4703226-01 

-.2 167026-01 

— .5105o0E-06 

-.1 lc^d46-06 

42 

.2452656*01 

.8702666-05 

-.2046826-05 

-.5052416-01 

-.1273106-01 

-.4657516-06 

-.7379716-08 

43 

. 3024666*01 

. 7678256-05 

-.164 7626-05 

-.5053196-01 

-. 3133o56-02 

3931626-06 

. 59 1 9&2E-0 7 

44 

• 309660E *01 

.7194506-05 

- . 1 02 7 80E -05 

-.4724o46-01 

•5C3u6 26-02 

-.3339156-06 

. 905516E-J7 

45 

. 316o69E*C l 

. 7140676-05 

-.7209076-06 

-.4125956-01 

. 131097E-01 

- . 2 G5Co5E-06 

• 1 236866-00 

46 

. 3240716*01 

. 7C6747E-05 

-.7759776-06 

-.3345026-01 

. 1 79366E-01 

-.22 249CE-06 

. 1 52723E-06 

47 

. 3312726*01 

.6046566-05 

-.931 1O7E-06 

-.2404036-01 

. 1 994 0 1 E-0 1 

-.1 5l606t-06 

. 1 59 75*» t-Oo 



48 .338474E+01 .6541006-05 -y 1072856-05 

49 .3456766+01 .6200726-05 -.117462E-05 

50 • 352677E+0 1 -585692E-05 -.124131E-05 

51 .36CQ79E+0 1 .5519SSE-05 -.128115E-05 

52 .367230E+01 ,5191056-05 -.12998BE-05 

53 .3744626+01 .4870156-05 -.130090E-05 

54 .3616846+01 '.454o52E-05 -.1284056-05 

55 • 3 C86 656 * C 1 .420277C-05 -.1241746-05 

56 .3 96087E + 01 ,3790668-05 -.114L45E-05 

57 .403288C+0 1 .317189E-05 -.7751846-06 

50 -410490E+01 .322463E-05 .105239E-05 

59 .4l7o91E+01 ,52469 7E -05 -.617916E-06 

60 • 424893E + 01 .442294E-05 -.113265E-05 

61 .4320956+01 .399100E-05 -.1200196-05 

62 . 43«29bE+Ql .3703726-05 -.120147E-05 

63 .4464966+01 .348C25E-05 -.1186266-05 

64 • 4 53o 99E +0 1 .329231E-05 -. 1165456-05 

65 • 46090 IE +01 .312813E-05 -.U4242C-05 

66 -4b8102E+0i .2580266-05 -.1113626-05 

67 .475304E+01 .2045O3E-O5 -,1094u4E-05 

68 .4825066+01 .2720136-05 -.1070G2E-05 

69 .489707E + 0 1 .2603956-05 -.1047366-05 

70 . 4969 09E + 0 1 .249552E-05 -.102^416-05 

71 . 504 1 10E + 01 .2393 73E- 05 -.1002078-05 

72 .511312E+01 .22976 0E-05 -.9004CoE-06 

73 .5185 L3E + 0 L .2206736-05 -.959471E-06 

74 .525 71 5E +01 .2 12030E-05 -.939320E-06 

75 .5329176+01 .203793Er05 -.92001QE-06 

76 .5401186+01 .1959116-05 -.90I629E-06 

77 .54 73 20E + 01 .188314E-05 -.8841566-06 

78 .5545216+01 . 1009406-05 -.867611E-06 

79 -5uI723E+0l .173724C-05 -.B5^055E-06 

80 .5639246+01 .1665326-05 -.8371116-06 

81 .5761266+01 • 1 5924 7E- 05 -.8223596-06 

82 • 5 C33 28E * 0 1 .151683E-05 -.0065026-06 

83 • 590529E + 0 1 .14361UE-05 -.786964E-06 

04 .5977316 + 01 .1346426.-05 -.7575G8E-06 

85 .6049326+01 .125559E-05 -.7060476-06 

86 .6 121 346 +0 1 . 1180106-05 -.616269E-06 

87 .619336E + 01 .U6361E-05 -.5109726-06 

88 .6265376+01 .1192006-05 -.445695E-06 

89 .633 7 39E + 0 1 .1Z2092E-05 -.4256086-06 

90 . 64 0940E +01 .123075E-05 -.429o6lE-06 

91 .640 142E + 0 1 .1222486-05 -.4393036-06 

92 .6553436+01 .120333E-05 -.4469736-06 

93 .66Z545E+01 .117878E-05 -.45115JE-06 

94 .6697476+01 .115108E-05 -.452225E-06 

95 .676948E+01 .1124266-05 -.450906E-06 

96 .68415CE+01 .109679E-05 -.4477976-06 

97 .69135 IE *01 .1070006-05 -.4433936-06 

98 .6985536+0 1 . 1044076-05 -.4381u8E-06 

99 • 705754E+C 1 .101912E-05 -.432J75E-06 

100 .7 12956E+01 .995228E-06 -.4262126-06 

101 .720152E+OL .972453E-Q6 -.419043E-06 

102 .7273596+01 .950042E-06 -.413405E-06 

103 .734561E+0L .930445E-06 -.4070306-06 


Lh 




• 164928 E— 0 1 • 1916236-01 

•930299E-02 .159458E-01 

• 397548 E-02 .109368E-01 

.917262E-03 .4941008-02 

• 220044E-03 -.1102236-02 

.1675266-02 -.663187 E-02 

.482474E-02 -.107602E-01 

.903813E-02 -.1314156-01 

.1360556-0,1 -• 1 36073 E-01 

.1703206-01 -• 1 2247BE-0 L 

■21 1 270E-0 l -.9380156-02 

.230653E-01 -.545J93E-02 

• 2343 36 F-0 1 -.116230E-02 

.2223956-01 .300973E-02 

.1969496-01 .6491806-02 

.1617436-01 .8877236-02 

.1215566-01 .992926 E-02 

.0151326-02 .959992E-02 

•464046E— 02 .8023656-02 

.200091 E-02 .548804E-02 

• 505494 E-03 .2386946-02 

• 2 18439 E-03 -.8370086-03 

. 107383 E-02 -k375J0ut-02 

.205500 E-02 -.5596296-02 

• 52371 OE-02 -.7312226-02 

• 7045 lo E-02 -.7o823OE-02 

• 1028u9E-0 l -.603172E-02 

• 122230E-01 -.5217776-02 

.1339916-01 -.3001906-02 

• 13.67576-01 -.50 96086-03 

• 130402 E-0 1 • 19 1452 E-02 

.1160146-01 .395463E-02 

. 956B70E-O2 .5361696-02 

.721959 E-02 .5982926-02 

•485925E-02 .577703E-02 

.2779 736-02 .481364E-02 

.122200E-02 .3250516-02 

.3469626-03 .134705E-Q2 

•2133046-C3 -.650055E-03 

• 798675E-03 -.2ho408E-02 

. L96942E-02 -.306343 E-02 

.350473 E-02 -.4682336-02 

• 519311 E-0 2 - • 48394oE-02 

.6779338-02 -.434470E-02 

.0042888-02 -.329364E-02 

.8816128-02 - • 1 E5105t-02 

.900454E-02 -• 226056E-03 

.859504 E-02 . 135740E-02 

.7656958-02 .2690776-02 

.6 32367 E-02 .3u0770F-02 

.4777036-02 .400420F-0 2 

.3222248-02 .3o49c7E-02 

. 105200E-02 .3187066-02 

.832900 E-03 .212795E-02 

.2731276-01 • 8 2 74056-03 

•210O2QE-O3 -.5300366-03 


-.8732086-07 .1430358-06 
-.3895496-07 . .109803E-06 
-.970809E-08 .689906S-07 

.126697E-0S .254490E-07 
- • 26751 9E -08 -. D &:>L 9 1 £-08 
-.1670O2E-O7 -.301 1 fabE-07 
-.3575446-07 -.4273C7c-07 
-.5430356-07 -.44 l07cc-07 
— • 671 C56E-07 -- 36O:jC8£-07 

— .66C573E-07 — • 2 502 5^ E-o7 
-•582551E-07 -.5246146-^7 
-.1244626-06 1 45oSroc- 07 
-.104WE-06 .2140146-07 
-.0514606-07 .3870356-07 
-.6514476-07 .477 Jb7t-j 7 
-.457o0bE-07 • 5COt>23E-J7 
-.2044008-07 .4605706-07 
-.1453136-07 • 39 J4206-07 
-.4834406-08 .291035c-07 

.292029C-09 .1701278-^7 

• 1 1 EG9EE-0S • 7034 096— u d 
- . 1 44 553E-0E -.195J956-O0 
-•O52436E-08 — • 8 2o563 c-06 
-.1264286-07 -.1149206-07 
-• l 52 C67E-07 1 Lo647E- 07 
-.245671E— 07 — • 92o* d9c-08 
-.282264E-07 — .432zo4c— j 6 
- .297 1 OCE -07 .61 19096-09 
- • 2895o8E -07 .bl9953E-J3 
-.26204CE-07 .1113146-07 
-.2 1 934CE-07 . 1477506-07 
-. lo7348E-0? .1673556-07 
-.114466E107 .1693516-07 
-.6576866-08 . I 54647c-,} 7 
- .271 C99E-08 .l2o322E-j7 
-.2030126-09 .91OJ41E-O0 

.8205716-09 • 53i953t-0d 

• 5 1 o5 1 56-09 • 193 od5E-08 

- .65 823U1E-Q9 -.632556E-j 9 
-.218 861E-08 -. 2459036-0 3 
-.4058736-08 -.3731^56-03 
-.6272226-08 -.4224SCe-„3 
- . e47O75E-O0 -.3724696-Od 
-.1015636-07 -.2333116-08 
-.1115046-07 -.3634026-09 
-• 1 1227JE-07 . 1755446-03 

104744E-07 .331 1056-08 
-•905lo96-08 .54C19oc-v»a 
-.71931 Cc-OE .6379956-08 
-.5166716—06 .66 o 41.2E-J8 
-.3233876-08 .627415E-08 
-.161 534E-08 .5316456-08 
-.4852546-09 .39o234c-0d 

. 8344d GE“ 10 .2419 u2 6-08 

.8235306-10 .899o43E-09 
-•419667E-09 -.4043476-09 



IQ* 

•741762E*01 

•91 13256-06 

-•4008566-06 

105 

.7489656*01 

.8935636-06 

-.3950436-06 

106 

• 7561656*0 1 

.8772646-06 

— .3 89798E-06 

107 

. 763367E*Q1 

.8627246-06 

-.3855426-06 

108 

. 770569E *0 1 

•P50028H-06 

-.3823086-06 

109 

. 7777 7QE *0 1 

.8392336-06 

-• 382462 E-06 

110 

• 7 849 72E+ 0 l 

.8301296-06 

3 05939E-06 

111 

•792l73E*0l 

.8214916-06 

-.393697E-06 

112 

• 799 3 75E * 0 l 

.8094636-06 

-.4142196-06 

1 13 

• 8 0c5 7oE* 0 1 

.7859886-06 

-.4411596-06 

114 

.8137786*01 

• 742o2 6£-0o 

-.4673216-06 

115 

.8209306*01 

.631 1906-06 

-.47501 7E-06 

116 

.8281816*01 

.6227616-06 

-• 45392 OE -06 

117 

i 835383E*Cl 

. 5 85o3 3E-06 

-.4 I94l26-0o 

118 

. 8425346*01 

.5654276-06 

-.3885286-06 

119 

.8497366*0 l 

• 5523 5 3C-06 

3o4052 £-06 

1 20 

• 6569386*0 1 

• 5423958-06 

-• 3‘*4o£05-0o 

121 

. 864 1896*01 

.5336476-06 

-• 329586E-06, 

122 

.8713916*01 

. 525289E-06 

-.3170136-06 

123 

.8785926*01 

.5170046-06 

- . 306351 E-06 

124 

.8057946*0 l 

.508695C-06 

-.2970576-06 

125 

.8929956*01 

.5003526-08 

-• 2 06 7696-06 

126 

.9001976*01 

.4919956-06 

-.2812396-06 

127 

. 9073996*01 

. 483 653E- 06 

-.2743016-06 

128 

.9146006*01 

.4753636-06 

— * 267 0o5E-06 

129 

• 92 1802t *C 1 

* 467 1 5 IE-06 

-.26181 4 E-06 

130 

.9290036*01 

.4590356-06 

-.2560766-06 

131 

.,’9362056*01 

.4514146-06 

-.251 1156-06 

132 

.943406t *0,1 

.4439846-06 

-.2464386-06 

133 

. 95060BE + 0 1 

.4365546-06 

-.2417626-06 

1 34 

• 957 d 106 ♦ 0 l 

.4291235-06 

-.2 3 7086E-06 

135 

.9650116*0 l 

.4216936-06 

-.232409E-06 

136 

.9722136*0 l 

*41 42636-06 

-.2277^36-06 

137 

.9794146*01 

.4068326-06 

-. 22305oC-06 

138 

.9366166*01 

.3994026-06 

-.2133806-06 

139 

. 99381 76*01 

.3919726-06 

-.213 7C4E-06 

140 

. 1001026*02 

• 38597 IE-06 

-.2101606-06 

1 4 1 

. 1008226*02 

. 3 801 74E - 06 

-.2068206-06 

142 

. 1015426*02 

.3743776-06 

-.2034616-06 

143 

. 1022626*02 

.3685806-06 

-.2001016-06 

144 

. 1029336*02 

.3627336-06 

-.1967426-06 

145 

.1037036*02 

. 3 56 90 76- Oti 

-.1933826-06 

146 

. 1044236*02 

•35119 CE -06 

-. 1900236-06 

147 

. 105143E*02 

• 345393E-06 

-.1866636-06 

148 

. 105863E *02 

.3395966-06 

-.1833046-06 

149 

. L06503E*02 

• 333799E-06 

-.1799446-06 

150 

. 1073036*02 

.3280026-06 

-.1705356-06 

151 

. 1080246*02 

•3222G6E-G6 

-.1732256-06 

152 

. 1087446*02 

.3164096-06 

-.16936LE-06 

153 

. 1094646*02 

.3106126-06 

-.1665076-06 

154 

. 1 10 18*t*02 

.3056026-06 

-.1636156-06 

155 

. 1 109046*02 

.3013496-06 

-. Itol 185E-06 

156 

. 1 11624F*02 

.2970966-06 

-.1587516-06 

157 

. 1 123456*02 

.2923436-06 

-.1563136-06 

158 

. 1130686*02 

.2885896-06 

-.1538046-06 

159 

. 1 137856*02 

• 2 04336E- 06 

-.1514506-06 


-•6467026-03 

-•176360 E-02 

-•129631E-08 

“• 134798E-08 

— • 14707 IE— 02 

-.2711746-02 

-.238543E-08 

-.18<>2UE-08 

-.2555906-02 

-• 326023 E-02 

-.3513036-08 

186380E-0B 

-.3738226-02 

-.3351706-02 

-.4517276-08 

-.19533iE-08 

—.4847 38 E-02 

-.2991726-02 

-.526566E-08 

-.687-.35E-09 

-.5729686-02 

- . 224568E-02 

- • 56b 742 6-98 

• 3 Co 7 ‘♦oc-09 

-.6 267 70 E-02 

-.1227766-02 

-.5676906-08 

• 1 40«jwioE-j 6 

-.6394536-02 

-.8399796-04 

-.5286296-08 

.24i>125t-0B 

-.6100696-02 

.1028356-02 

-.4512326-08 

.3359436-08 

-.5433246-02 

. 1562C6E-02 

-.3404086-08 

.3939086-^8 

-.4487466-02 

.2599136-02 

-.2117866-08 

• 40272?c-u 0 

-.3392396-02 

.2 S6536G-02 

-.9497906-09 

• 3 5o332E-US 

-.2292446-02 

.2733656-02 

1 84 5 1 4 1-09 

• 2740UE-J8 

-.1327706-02 

•225000E-02 

.1661346-09 

• 187453 E-G 6 

-.6153916-03 

.1477536-02 

.2261036-09 

.1074536-08 

-.2351546-03 

.5344206-03 

• 6466d(»C- 10 

.3607976-09 

-.2 197396-03 

-.4474006-03 

-.2734856-09 

— • 1668o*E— u9 

-.5522956-03 

-.1335426-02 

-.7348676-09 

-. 5306166-09 

-.1 170326-02 

-.2014226-02 

-. 1253&9E-08 

-.6870*26-09 

1975546-02 

-.2400236-02 

-.1756*86-08 

-.6357216-09 

-.28482OE-02 

-.2451776-02 

- • 2 L 771 P.E-06 

-.4011266-09 

-.3663736-02 

-.2173196-02 

-,24607Uc-08 

-.2933656-10 

— .430950 E-0 2 

-.1612736-02 

-.2573816- 08 

• 4 1 6:>4oE-09 

-.4699396-02 

-.8546656-03 

-.2507316-08 

• B 75t>d 5E-09 

-.47344 2E-U2 

68^1796-05 

-.22761 76-08 

. 1278326-08 

-.455802 F-02 

.8143896-03 

- • 1 91 807E -06 

.1573796-08 

-.4055556-02 

.1500056-32 

- . 1 4775 IE-00 

. 1727 1 Ofc-O 8 

—.334821 6-02 

.1962756-02 

-.1018556-08 

. 1 72od0t-08 

-.25324 IE-02 

.2147686-02 

-.5950776-09 

. 1577626-08 

-.1716096-02 

.2038676-02 

- .2562466-09 

• 1 3C*9o t-u B 

-. 100400 6-02 

.1660116-02 

-.37251oE-10 

• 9 50*246-09 

-.403749 b-03 

. 1070506-02 

.448C01E-IC 

• 5 635*56-0 9 

-.2147196-03 

•3 95506 E-03 

-.7989776-11 

. 19o l 7 l c- 09 

-.2209666-03 

-.335287E-03 

175830E-09 

-. 107*5V6-u9 

-.409 L 116-03 

-.105189E-02 

-.4250656-09 

-• 3 133 iuE-09 

-.9712476-03 

-.15572 9E-02 

-.7134996-09 

-. 4C2 a 53 1-09 

-.1592396-02 

-.1838626-02 

-.1001066-08 

— .3 749ooc-09 

-.226 131 E-02 

1864916-02 

124639E?08 

-.2*1 3076-09 

-.2883166-02 

-. 1640206-02 

-. 1413Ue40B 

-.27**536-10 

—.337231 E-02 

-.1202016-02 

-. 14S349E A 90 

. 2 31 7uo t-09 

-.3663456-02 

-.6152946-03 

- • 1 45C09E -08 

• 4 975 3o 6-J9 

-.3719796-02 

.3698776-04 

-. 1 3 2 C 7 6 E - C 0 

• 7 323*8 56-u9 

-.3537076-02 

. 6o55 1 6E-C3 

-.1115726-08 

.9056*66-09 

-.3 14324E-02 

. 1186736-02 

-.8641346-09 

• 99o5 14C-09 

-.2593796-02 

.1533616-02 

-.5996876-09 

.9963266-09 

-.1963526-02 

•1 o65176-02 

-.355784E-09 

.9051 56E-09 

13 3o02 6-02 

. 1569276-02 

-.1611096-09 

• 75o^*tc-09 

-.7922256-03 

. 1265346-02 

- • 3t>06v 76- 1 C 

• 54*9j*C-^9 

-.3996526-03 

.7 999 3 5 E-03 

.942 847E-1 1 

.320 '>9*6-09 

203760 E-03 

.2399356-03 

-.2333566-10 

. lCd*546-09 

— .222466 E-03 

- ,33o943t-03 

-.1231166-09 

-.6637096-10 

— .4*450 26-03 

-.8523516-03 

-.2 714L7E-09 

IG33^9E-o9 

-.8317226-03 

-. 124023E-02 

- .443990E-09 

-• 23at31 E-U9 

-.1324936-02 

-.1450596-02 

-.0147466-09 

-.21 7oS*6-09 

-.1052386-02 

-.1460446 -02 

-.7593156-09 

L3o*L8E-09 

-.2339706-02 

1273836-02 

-.8591 826-09 

78490*6-1 l 


160 

•114505E+02 

.2800836-06 

-.1490168-06 

161 

. 1 152256+02 

.2 75830E-06 

-.1465B2E-06 

162 

.1159456+02 

.2715766-06 

144148E-06 

163 

.1166668+02 

.2673236-06 

-**.141 7148-06 

164 

.1173 Q6E + 02 

• 263070E-06 

-.1392818-06 

165 

. 1 18106E+02 

.2508176-06 

-.1368478-06 

166 

. 1 18026C + 02 

. 2 54563E- 06 

-. 134'»1 3E-06 

167 

. 1195468 + 02 

.2505878-06 

-.1321368-06 

1 08 

. 120266C + 02 

.2474436-06 

1303316-06 

169 

.1209068+02 

•2442996-06 

-.1285258-06 

170 

• 12 17076* 02 

.241 1556-06 

-.1207198-06 

171 

. 122427E + 02 

• 238G 1 16-06 

124914E-0o 

172 

. 1231476 + 02 
i 123967E + 02 

• 234 866E-06 

-.1231086 -06 

173 

.231 7226-06 

-.1213036-06 

174 

. 1245076 + 02 

.2285706-06 

-.1 19497C-06 

175 

. 125307C + 02 

•2254346-C6 

-.1170926-06 

176 

. 12o023E + 02 

.2222906-06 

-. 1153868-06 

177 

. 126 740E + 02 

.2191456-06 

1 140J1E-06 

178 

. 1274686+02 

.2160016-06 

-.1 122756-06 

179 

. 12G188E + 02 

.2123576-06 

-.1 104706-06 

100 

.1289036+02 

.2097136-06 

-. 1 03664E-06 

181 

. 1296288+02 

. 2 0656 86 — 06 

IO605OC-O6 

102 

. 1 3034°E + 02 

.203424E-06 

-.1050536-06 

183 

. 1310696*02 

.2003338-06 

-.1035916-06 

184 

.1317096+02 

.1985136-06 

1022266-06 

185 

. 1325096+02 

. 196142E-06 

-. l00bblE-06 

136 

.1 33229E+02 

.1937726-06 

-.9549606-07 

187 

<1 J3949E+02 

.1914018-06 

-.9813 ICE-07 

168 

. 1346696+02 

.1890318-06 

-.9670616-07 

199 

.1353908+02 

.1866606-06 

-.9540116-07 

190 

. 13ol 108+02 

• 1 04 29 C6- 06 

-.9403618-07 

191 

. 1368308*02 

.1819198-06 

-.9267118-07 

192 

.1375506+02 

.1705496-06 

-.9130618-07 

193 

.1302706+02 

.1771788-06 

-.8994126-07 

194 

. 1309 90l ♦ 02 

.1740086-06 

-.3857628-07 

195 

• L3971 16 + 02 

.1724376-06 

-.0721126-07 

190 

. 1404318^02 

.1703306-06 

-.0599866-07 

197 

•14115 1E+02 

• 1 604o 7E-06 

-.8492058-07 

198 

•14137 16 *02 

• 16660 3E - 06 

-.0335456-07 

199 

. 1 4259 1 E ♦ 0 2 

.1647396-06 

-.8273258-07 

200 

. 14331 16*02 

.1628766-06 

-.8171056-07 

201 

. 144032E+02 

.1610128-06 

-.80o3Q4E-07 

202 

. 1447526+02 

.1591496-06 

-. 795O04E-07 

20 J 

. 1454728*02 

.1572058-06 

-. 784944E -07 

204 

. 1461926*02 

. 15542 IE- 06 

7742246-07 

205 

. 1469128+02 

. 1 53 5586-06 

7635C3E-07 

206 

. 1476 328 *02 

.151694E-06 

-. 752 78 36-07 

207 

.1483528*02 

. 14963 IE-06 

7 42 Go 3 E-0 7 

208 

.1490736*02 

.1479676-06 

-. 731 3436-07 

209 

. 1457936+C2 

• 1 4o20 16-06 

-.721 1926-07 

210 

.1505138*02 

.1447128-06 

-.7120048-07 

2 1 1 

. 1 512336*02 

.1432226-06 

-. 7C4l3oE-07 

212 

. 1 519538*02 

. 141 7336-06 

-.6556076 -07 

213 

. !52o73t *02 

.1402448-06 

-.6670796-07 

214 

.1533948*02 

.1387556-06 

-.6785518-07 

2 15 

. 154 1 148*02 

. 1 372668-06 

-.6700238-07 


-0 

-4 


-.271994E-02 -.9210066-03 

-.29423 l E-02 -.453713E-03 

-.2978616-02 .6232276-04 

-.2026328-02 .5565598-03 

-.2508378-02 .963213E-03 

20t>933 E-02 .122995E-02 

-. I 36892 E-02 .1324558-02 

-.1073608-02 .123870E-02 

-.6475 158-03 .9bR5l06-03 

— •’343914 E— 0 3 .611/386-03 

-• 198324 6—0 3 • 1 62 1 62 F-03 

-.224218 fc-0 3 -.2900131-33 

-.41 1 760C-0 3 -.7067118-03 

— . 729699 E-03 -.101030E-02 

-.1 1 3009 E-02 -• 1 1 70 1/9 E-02 

-.1555146-02 - . 1 169558-02 

1945196-02 -.101 1266-02 

-.2246 75 E-02 - • 7 20 7 80E - 03 

-.2419:>b6-02 - . 3 403 *.9 E-03 

-.2441026-02 ■ 7 t»7523G-04 

-.2 3 1 1 77 E-0 2 .4735558-03 

-. 2049266-02 *7973236-03 

— . 1690 788-02 . l'C06l5£-02 

-.1295026-0 2 • 107473 E-02 

-.885946E-03 .997105E-03 

-.5453708-03 .7871396-03 

-• 306097 E-03 .476129E-03 

-.1963146-03 .1002248-03 

-.22622 1E-0 3 -. 2 c5 799E-03 

-.38 70 30 C-0 3 -.5'95505E-03 

-.6526938-03 -.8375816-03 

-.9033056-03 901453 E-03 

-. 1331 73E-02 -.9 531378-03 

-• 1049 14 E-02 -.8 10 7288-03 

-.1892108-02 -.573444E-03 

-.2028 17 E-02 -.258392E-03 

-.204003E-02 .844571E-04 

-• 192745 E-02 .4U6309E-03 

— . 170696 E-02 -O701 37E-03 

1409? 7 E-02 .830J55E-03 

-.1074308-02 . 8 0 57 53 C-03 

-.748049T-03 .G1512GE-03 

—.471643 6-03 .6 3o2 8 5 E-03 

-. 2 80 1 49 fc-0 3 .3 750296-03 

— • 190552 E-03 .699932C-04 

-• 22 84 7o £-0 3 - . 2 335 33C-03 

-.3679608-03 -.5033708-03 

-.5927o3E-03 - . 70402 Ot- 03 

— .8696/4 E-0 3 -.800 7478-03 

-.1159178-02 - « 78 7 599E -03 

-. 142039 1-0 2 -.6686658-03 

-.lo 1909 E-02 —.4621 64 E-0 3 

-.1727326-02 -.1977248-03 

— .173 IS0E-02 . 0 789 796-04 

-. lo 3 30 t E-02 . 3557b2E-03 

-• 1 445 34 L -02 .5701 1 G E -03 


-•8990526-09 .1473566-09 

878C846-09 .3061436-09 

-.7999356-09 .4462666-09 

-.6766686-09 .54531 U-G9 

-. 525 7 19E-09 .602799E-G9 

-.3672628-09 .6C15l2E-^9 

— • 221 363t-09 • 54b004E-Q9 

-.1053598-09 .4522ojfc-J9 

-.3139038-10 .3269916-09 

- • 539422E-1 1 . 1 93o4 1 1-09 

-.2 72 777E-10 .64237o6-LU 

-.9059258-10 -.4292238-10 

-.1C37UE-09 -. 1 152926-09 

-.291t406-09 -.1455*68-09 

-.3981 97E -09 -.1325298-09 

- .4GC22CE-09 EC«2 556-l C 

-. 5495058-09 .62di>2G6-U 

- . 574 592 E -09 *983545l-10 

- . 5oGd3 58- 09 . 19614CE-09 

-.5112368-09 .28fc^u28-09 

-.4333446-09 • 3 5u:>2 5t-09 

-.33511 38-09 .383oOJ8-09 

- • 2 3324 5E-U9 .3 922^76-09 

-.14680/1-09 .34401 l c - j9 

-. 73941 88-10 • 2 6 6 5 C 5 1 - 0 9 

-.275 8016-10 .203956-09 

-.1194006-10 .1227166-09 

-.2695476-10 * 3597o£,t-lQ 

-.68463GE-10 -.283537E-10 

-.1290558-09 -.7423466-10 

-.1990298-09 9 295^36- 1 0 

- • 26 79t> l 6-09 -.8373276-lC 

-. 3261 3 7E-0*; -.495395E-10 

-.3656458-09 .3619018-11 

-.3615478-09 *&73524E-iO 

-.3722666-09 .13232J C -G9 

-.3402158-09 .189825E-J9 

- . 290C 2 58t C9 *2324Sob-J9 

-.22dld7Ef09 *2547928 -J 9 

162552t L 09 .25439-.E-U9 

-.1026848-09 • 2 32 JvCi- j9 

-.54714 58- 10 . 191 ^o7G-o9 

-.244 1866-10 . I 3b 763 l-J9 

-. 1457046-10 « 6 lu7 1 0E- 1 J 

- .2 51 293E-1 0 .2604598-10 

-.53256 38-10 -. 191d 441:-1 j 

-.9408678-10 -.4541 7^8-1 0 

-.1410566-09 -.6 It >748-10 

-. 1872456-09 -•9*»b3l3L-10 

-.2262 726-09 - . 3 1 5-* 94 6 - 1 0 

-.2532738-09 .44,Gj5l-11 

-.2644338-09 .473136 t- 10 

-.2585728-09 .9212448-10 

-.2368388-09 .1313io8-09 

-.2024528-J9 . I Gul 63 E-G9 

-.160197E-C9 .1750548-09 



216 

• 154834E+-02 

•135777E-06 

-.661494E-07 

217 

. 1 55554E+-02 

• 1 34287E-06 

-• 652966E-07 

210 

. 1 562746+02 

.1327986-06 

-. 644438E-07 

219 

. 1569946*02 

.1313096-06 

-.6359096-07 

220 

. 1577156 *02 

. 1298206-06 

-.6273816-07 

221 

. 1584356+ 02 

• 1 2 8 33 IE-06 

-.6188536-07 

2 22 

. 1591556+02 

.12684 IE-06 

— • 6 1032 5E-0 7 


******** WARNING. 

EXTRAPOLATION REQUIRED FOR * I TJ S. 

i 


-.1194856-02 -703343E-03 

-.915521 E-03 .7391066-03 

—.64455 1 E— 03 .6745856-03 

—.4 L7258 E-03 .5205086-03 

2 (>2380 £-0 3 .2994376-03 

-.1983446-03 .4249306-04 

-.230990E-03 -.215005E-03 


-.11 5707E-09 .1745366-09 
-.7466186-10 .1590336-09 
-.4212246-10 .131121 t-09 
-.2169016-10 .9480126-10 
-.1527606-10 .5533416-10 
-.2282396-10 .1772776-10 
-.424214E-10 -.1317376-10 





J 

TIME* TJCJI 

1 

0 , 

2 

• 3000000E-01 

3 

•6000000E-0! 

4 

•90000006-01 

5 

. 12000C0E+00 

6 

. 15000006*00 

7 

. iaoocco 6 *oo 

0 

• 21000C0E* 00 

9 

• 24 J0C00E*Q0 

10 

• 2 7000CCE*00 

11 

• 3C30QCCE *00 

12 

• 33000CGE*QQ 

13 

• 3OG0QCGE*00 

14 

. 390Q0 COE* 00 

15 

.42000006*00 

16 

• 45000CCE *00 

17 

•4R0O0CGE*00 

18 

. 51000006*00 

19 

•54300006*00 

20 

. 5700CCCE*00 

21 

.60000C0£*00 

22 

.63000006*00 

23 

• 60 OOOOOE+OO 

24 

. 6900000E+00 

25 

- 72000C0E+00 

26 

•75000006*00 

27 

. 78000C0E * 00 

28 

•81000006*00 

29 

.84000006*00 

30 

. 87000C06*00 

31 

•9900000E+00 

32 

• 9300000E+00 

33 

. 9600C COF * CO 

34 

•9900000E*00 

35 

. 102 OOCOE*01 

36 

.10500006*01 

37 

• iocooqoE +01 

30 

.11100006*01 

39 

.11400006*01 

40 

. 1170000E*01 

41 

. 12000C0E+C1 

42 

- 12300006*01 

43 

. 1260000E + 01 

44 

• 1290000E + 01 

45 

. 13200C06* Cl 

46 

. 13500CCE + 01 

47 

. 13800006*01 

48 

. 14100006*01 

49 

• 14400006*01 

50 

. 14700006*01 

51 

. 15000006*01 

52 

. 15300006*01 

53 

. 15o00006 *01 


RESPONSE FUNCTION, X(J) 


• 6993698E-07 
-.26589856-07 

— *49 72459E-06 
-• 1793572E-05 
-■4321 107E-05 
-■ 8406048E-05 

1424567E-04 
-.2180131 E-04 
-.3094044E-04 
-.4141 172E-04 
-• 5208136E-04 
— • 64 50292 E-04 
-• 7b992 08E-04 
-.8866676P-04 
—•9945 10BE-04 
1088571E-03 
11656U5E-03 
-. 1223 047E-03 
12 5918 76-03 
12 73199E-03 
126491 lb-03 
1236734E-03 
-.1191 339E -03 
-. U32369E-03 
106422 7fc-03 

— .9922 701E-04 
-.92213016-04 
-. 05920 75E-04 
-- 8083034, E -04 
-.77199546-04 
-. 7487070E-04 

735652 1 E-04 
72963196-04 
-•72720J2E-04 
724638 7 E-04 
71 82 909E-04 
-.70405326-04 
-.60134076-04 
-.64505936-04 
-.59390586-04 
-. 52745606-04 
-.445 8186 E-04 
-.34992456-04 
-.241141 3E-04 
12204836-C4 
. 4202 769E -06 
. 1345 L 1 2E-04 
. 265L3806-04 
.3921719b *04 
.51 1 9703 E-04 
.62127936-04 
.71 73565E-04 
. 7979629E-04 


1 

TIME, TUI) 

FORCING FUNCTION, Fit I 

l 

0. 

0. 

2 

. 135604 0E-01 

. 10371866*00 

3 

•2712080E-01 

•20631 85E*C0 

4 

.40681206-01 

.30669296*00 

5 

.54241606-01 

.40375906*00 

6 

• 6 7802 OCE-Ol 

-49647006*00 

7 

• 8 13624 CL- 0 L 

.58392586*00 

8 

.9492280C-CI 

• 6648 E 406 *00 

9 

. 10d43J2c*00 

.73877036*00 

10 

.12204366*00 

• 80,60 ?dc *00 

11 

. 135604CE*00 

• C61 92 55c *C0 

12 

. 14916446*00 

. 9058o5 Jc* OC 

13 

. 162 724 86 * 00 

•947991 7t*00 

14 

.17628526*00 

• 975392ot *00 

15 

. 1898456t*00 

.99326576*00 

16 

. 2034 06 OE *00 

• 99992 54c* 00 

17 

• 2 169604E *00 

• 995 799.>b + 00 

10 

. 23C520 b6* 00 

• 98093 1 at* GO 

19 

.24408726*00 

• 9 554 034c ♦ CO 

20 

• 25764 76E *00 

.91972356*00 

21 

.27 1 2 00 CE *00 

.87405286*00 

22 

• 2 04763 46*00 

.3,1894906*00 

23 

• 298323Bb*00 

.755011^6*00 

2 4 

.31 16692E + U0 

.68293056*00 

25 

.32544960 00 

• 603 4 b2u£ * 00 

26 

.33901006*00 

.51 7524*6*00 

27 

.33257046*00 

• 425954d£ *C0 

28 

.3o6l3086*00 

.329649^6*00 

29 

• 379o9 1 2t *00 

. 2 30 1 567u*w3 

30 

.39325166*00 

•1279813c *00 

31 

. 4 068 12 OE *00 

• 2442 33lt— 0 1 

32 

.42037246*00 

- • 79394 12E-C1 

33 

• 4 33932 8E *00 

-• 16235 7it*00 

34 

• 44749 J2t *00 

-.2 33 3531 E *00 

35 

• 4o 1 033 tE*00 

-.33129276*00 

36 

• 4 746 1 4CF*00 

-• 475 1 1 S*»t * 00 

37 

. 4 u 8 1 7,46*00 
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38 
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-.64o44i2b*00 

39 
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40 
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41 
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42 
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43 
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44 
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- . 9 7C2 7 5^6* GO 

45 
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-•990 139>c*U0 

46 
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- .99932886*00 

47 

.62377046*00 

- . 997 73 SoE ♦ 00 

48 

• 6373 33 OE * 00 
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49 
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-. 962 4 0456 *00 

50 

• 6644396b *00 

-•929 34236 ♦ CO 

51 
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52 

.69158046*00 
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53 
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• 6063097E-04 
• 5970552E-04 
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.22374666+01 0* 

.22510266+01 0. 

.22645876+01 0. 

.22781476+01 0. 

.22917006+01 0. 

.23052686+01 0. 

.23188286+01 0. 

.23323896+01 0. 

.23459496+01 0. 

.23595106+01 0. 

.23730706+01 0. 

.23866306+01 0. 

.24001916+01 0. 

.24137516+01 0. 

.24273126+01 0. 

•2440872E+01 0. 

.24544326+01 0. 

• 2467993E+01 0. 

•24815536+01 0. 

.24951146+01 0. 

.25086746+01 0. 

.25222346+01 0. 

•253579 5E+01 0. 

.25493556+01 0. 

.25629166+01 0. 

.25764766+01 0. 

.25900366+01 0. 

.26035976+01 0. 

.26171576+01 0. 

.26307186+01 0. 

.2644278E+01 0. 

. 2657838E+01 0. 

.2671399E+01 0. 

.26849596+01 0. 

.26985206+01 0. 

.27120806+01 0. 

.27256406+01 0. 

.27392016+01 0. 

.27527616+01 0. 

.27663226+01 0. 

.27798826+01 0. 

.27934426+01 0. 

.28070036+01 0. 

.28205636+01 0. 

.28341246+01 0. 

.28476846+01 0. 

.28612446+01 0. 

.28748056+01 0. 

.28883656+01 0. 

.29019266+01 0. 

.291 5486E+0 l 0. 

.29290466+01 0. 

.29426076+01 0. 

• 295616 7F +0 l 0. 

• 296972 8E +0 1 0. 

. 298J280E +01 0. 
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• CONVOLUTION NUMBER 2 • 



I 

FREOUENCY 

INTERPOLATED 

RESPONSE 

INPUT FREUUENCY 
FUNCTION 
FRF 

TRANSFORM OF 
IMPULSE TIME 
T FRF 

THE GUST 
FUNCTION 

TRANSFORM OF THE 
TIME CONVOLUTION 
FRF • TFRF 



REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL IMAGINARY 


1 

0 . 

-.694058E-03 

• 1 53553E-03 

-.175538 E-03 

0 . 

• 121 834E-06 

-•269545E-C7 

2 

• 7 20158E-01 

.2377216-01 

-.1505006*00 

.8740436-02 

.4750666-01 

•73575L6-02 

-.18ol00E-03 

3 

• 1440326*00 

♦364C49E-01 

-.4282516*00 

.3464026-01 

.8948676-01 

.3956386-01 

-. 1157 70c- u 1 

4 

* • 2 ’ 60476*00 

-.2521136*00 

-.8294676*00 

•750579E-0 1 

.1209216*00 

.813 771E-G1 

-.927430E-O1 

5 

•203O&3E+OO 

-•835076E*00 

-.8160716*00 

.1261376+00 

• 1 3 7759E+00 

* .7086736-02 

-• 2 17-*7o£+00 

6 

■ 360079E+00 

-.1 10383E* 01 

-.4965986*00 

• 1 82936E +00 

.1372046+00 

-.1337936*00 

-•2‘*24*7c + 00 

7 

. 4320956*00 

-• 1 17394E*01 

-.2543536*00 

.2401276+00 

.1 10JO2C+OO 

— .2 51 7566 * CG 

. 2Gv> 02 78 + 00 

e 

.5041106*00 

-.1204506*01 

-.9792886-01 

.2919926+00 

.8153 78 E-0 1 

-.343,7206*00 

I2e>6'07t+ JO 

9 

.5761266*00 

-.1243766*01 

.240711E-01 

.3334296*00 

.209740E-J1 

-.415^046*00 

-.2 801 09c“J l 

10 

•648l42t*00 

-.1306P7E+01 

. 135712E*U0 

.3501606*00 

-.3 574886-01 

-. 4o5 6306 + 00 

.9559716-01 

11 

.72Ci586*0C 

-.1397136*01 

.2507106*00 

.3691536*00 

— • 10795OL+00 

-.4870206+00 

« 2 463496 + 00 

12 

. 7921736*00 

-.1525436*01 

.4175526*00 

•35387uE+00 

-.1823536+00 

-.4713006+00 

• 42t01 7c+uC 

13 

.8641396*00 

-.1702276*01 

• 65604 1 E*00 

.3294156*00 

-.'2534006*00 

-.3945116+00 

• 64 74056* JO 

14 

.9362056*00 

— .1 92 5S7E* 01 

• 105 765E + G 1 

.282**31t*00 

-•315659c +00 

-.2099166*00 

.907094E+JJ 

15 

. 1008 226*01 

-.2116696*01 

• 1 60J89E+0 1 

• 220992 E+00 

— • 3o5384 E *00 

•19316 E 6*00 

• 1 173 lob* J 1 

16 

.1030246*01 

-.1794036*01 

. 3 lob loO+O L 

.1492706*00 

- . 3405OGE *00 

. 99392 1 1 +0C 

. 1 1 S753c+G1 

17 

. 1152256 + 01 

.29t>667E*00 

.4252916*01 

• 72 1545E-0 1 

— • 4 l 32 53 E + 00 

.1778966*01 

•1 S42o7c+00 

18 

’. 1224276*01 

m2 14826E*01 

•269391 E*01 

-.5188586-02 

-.4091596*00 

.1091096+01 

-.8929596+00 

19 

'• 1296286*0! 

• 2 095 1 5E *0 1 

.1 19389E*01 

-.7776336-01 

-.3671 b2E*00 

.2993256400 

-• 90hJ4**c + JO 

20 

..1368 306*01 

.1633016*01 

.5000786*00 

-.1411016*00 

-.'3496366*00 

— • 52 874 7E-0 1 

— • o44t» 79E+ 0 J 

21 

. 1440 J2E *0 l 

.1254336*01 

.2143996+00 

-• 19 1903 E*00 

-.2999056*00 

- . 1 764 1 Oc + 00 

-.417 3236+00 

22 

.1512336*01 

.9832246*00 

.7812466-01 

-.2276896*00 

-.2421006*00 

-•204955E *00 

-.2558346* 30 

23 

-15843pE*01 

• 79023 IE *00 

•11 6980E-01 

-.2475216*00 

-.160/026*00 

-. 193*»35E +0C 

-• I45b52fc+CG 

24 

.1656366*01 

.6454396*00 

-.2348896-01 

-.251709E*00 

-.120C8oE+00 

-.165263E+00 

-.7ip973c-Jl 

25 

. 1 728 38E *01 

.53324 IE +00 

-.4160236-01 

-. 24 l 75 l E*00 

-.642398C-01 

-.1315696*00 

-.241 7obc-Jl 

26 

• 1 00039E *01 

• 443540E *00 

- • 4978o0c-0 1 

-.2201246*00 

-.1539136-01 

-.98450CE-01 

• 3a8b9 1 t-ui 

27 

• 1 87 24 IE *0 1 

.3695336*00 

-.511496E-01 

-.1899836*1)0 

• 21 io 19 E-0 1 

-.691 234E-01 

. 175t52E-Jl 

28 

. 1944436*01 

.3067266*00 

-.4685416-01 

-. 1548546*00 

•4720116-01 

-.4528256-01 

• 2 1 75 796-01 

29 

• 20l644E*0 1 

• 25 1 25 IE * GC 

-.3650796-01 

— .11 0 245 E +00 

.61 7944E-01 

— .27453 It— Cl 

. 1 9o4^ 8 1-31 

30 

• 20E046E * 0 1 

.2009726*00 

-• l 74 7 5 CE-01 

-.8337o4C-01 

.656041E-01 

-.156049E-01 

. 1 4o4 1 66-0 1 

31 

.2160476*01 

.1572046*00 

.1647616-01 

—.529006 E-0 1 

.6C4085E-01 

-.931 I46E-G2 

• 8624658-^2 

32 

• 2232496 *01 

.1323426+00 

.74725JE-01 

-• 2b 7 Lo 5 6-U 1 

• 4 L 43 12 E-0 1 

— . 7*t2 3 1 6E-02 

• 4 2 7026 t-o2 

33 

• 2 3045GE + 0 1 
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• 1 47 3 3 7E +00 

-.1187456-01 

• 323957E-0 l 

— • O75947E-02 

• 36o960c-j2 

34 

.2376526*01 

.2601406+00 

.143L70E+C0 

-.257132E-02 

.1474676-01 

2 7dCl 86-02 

•34o8C6c-J2 

35 

•2443546*01 
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-.227595E-03 
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36 
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•L02O746-01 

-.3638066-02 

U24I 7E-01 

-.844745E-03 
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37 
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.2464086 *00 

-.2 153026-0 1 

-.III 080 E-OI 
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38 
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-.209766 E-0 1 

-.3167066-01 
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t C7C 6u£-u2 

39 
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-.4703056-01 
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-.3225706-01 

-. 74933 7E-02 

-.47:097E-J2 

40 

. 2 8066 1 E * 0 I 

.1684206*00 

-.522 9486-0 1 

— • 40 342 7 E-0 1 

-.2O5969E-01 

-.8269976-02 

-. 2 7Go 5 7 c- 02 

41 

• 28H063E *01 

.1467636*00 

-.5405146-01 

-.470322 E-0 1 

-.2167026-01 

-.809122E-02 

-• oC0594E-03 

42 

-295265E+01 

.12 51 5ut + 00 

-.5371956-01 

-.5052416-01 

-.1273 10E-01 

700727E-02 

.1 12 J78fc-02 

43 

• 302*1 666 + 01 

• 1 C4 9 6 IE * 0 0 

-. 45924 7E-01 

-.5053 L9E-0 t 

-.3I33o5E-02 

54**7 76C-G2 

. 1991 736-02 

44 

• 309o6L>L* 0 1 

.9387826-01 

-.317544E-01 

-.4 72464 E-0 1 

. 5 J3o62t -02 

-.4250076-02 

• 2049 2 2 c - j 2 

45 

• 3 16069F + 0 1 

• 92 043 IE-01 

-.2 1G5 1 76-01 

-.4125956-01 

• l j 1097 E-0 1 

-.354419E-C2 

.2113736-02 

46 

.324071E+C1 

.9347406-01 

-.1915556-01 

-.334532 F-0 1 

.1 79366E-01 

- • 2 783 14E-C2 

.231 7368-02 

47 

.3312726*01 

.9198946-01 
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-.2404836-01 

.1994616-01 

1B9G46E-02 
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- • 5 3boo3c-0o 

-.4117606-03 

-.7067116-03 

-.2 440j3E— 05 

-.1 50o42£-05 

7296996-03 

10 1030E-02 

•-.3e7o:cs-05 

-.1 89o49c-j3 

1 13U091-02 

-.1170596-02 

.-.5277016-05 

-.1 713376-05 

-. 1555 146-02 

-.1169556-02 

64615 7E- 05 

- . 101 bo4c-o5 

1945196-02 

-. 101 1266-02 

72o4 l 1 E-05 

.611 732E-J7 

-.2246756-02 

-.720 7 80 1-03 

-. 75 *566 - 35 

« 1 334086-05 

-.2419366-02 

-.3403496-03 

- - 7 392 u4 c-05 

.2671336-03 

-.2441626-02 

•7u7328F-04 

-.6728246-05 

. 3823366-05 

-.2311776-02 

.4735936-03 

-. 5b9247c-05 

• 4 672072-^3 

-.2049266-02 

.7973236-03 

-.4430056-05 

.5100*76-05 

— . 1690 78 E-0 2 

. 1006156-02 

-.3111516-05 

.5072126-05 

-. 123502E-02 

. 1074736-02 

1 9CG 7cE-05 

• 4623 54c- J3 

-.6839*oE-03 

.9971036-03 

-.9495746-06 

•3317356-05 

-.5453736-03 

.7671*96-03 

-.3445966-06 

• 276bo7c-w»5 

-• 30o09 7 6-03 

.4761296-03 

-.1464456-06 

. 1621 71 c-05 

-.196314E-03 

.10822*6-03 

-.3523026-06 

• 53u354c-0o 

— .22622 1 003 

-*2657996-03 

- • 906 8636-06 

-• 3681 C2 t-Oo 

-.3870306-03 

-.395309003 

-.171 G72E-05 

-.96651 3 c-»/6 

-.6525936-03 

-.C37pbit-03 

-.2636116-05 

-• 1205c36.-w3 

-.9833056-03 

-.9614536-03 

-. 33433c6t05 

-.1002116-05 

1331 73E-02 

-.9331376-03 

-.4309756-05 

-• 02492 9c-0o 

-. 1649146-02 

-.0167266-03 

-.482626E-05 

• 80 1 6 v2 6-0 7 

189210 E-0 2 

-.5734446-03 

-.5030126-05 

.5226026-06 

-.2028 L76-02 

-.2503926-03 

-.4501556-05 

.1779016-05 

-.204003002 

.8*43716-04 

-.4i73356-05 

•253*636-^5 

1927456-02 

.4033696-03 

- . 36C 72^6-05 

• 3C532 56-03 

-.1706966-02 

.6701876-03 

-.2559116-05 

.3381596-35 

1409276-02 

.0360556-03 

-.2128306 -05 

• 3370686-35 

1074806-02 

.8857536-03 

-.1334976-05 

.3071 156-05 

-.7400496-03 

.8151266-03 

-.7055416-06 

. 2 5322 3 6-0 3 

-.4715436-03 

.6362056-03 

-.310355E-06 

. lb*31Cc-J5 

— .2G0149E-03 

.3756296-03 

-.1854646-06 

.1070756-05 

-.1963526-03 

.6599326-04 

-.3293006-06 

• 34o3 7E E-*jo 

-.2284766-03 

-.2305336-03 

- • 7C4c60t-96 

-- 2* 65 96 c 6 

— . 36 79uOE-C 3 

-. 5U63 706-03 

-.1245156-03 

- . 643-# * * £ -0 6 

-.5927o3 E-03 

-.7040206-03 

-.1863646-05 

799*37 £ - 0 o 

-. 669o7* E-03 

-.8007476-03 

-.2 474136-05 

-• 70o3 7 L 6 - J6 

-.1139176-02 

-.7573996-03 

- . 2586 766-05 

— .3974 ^3c-0o 

-.1*20096-02 

— .6o8oc3fc-03 

-.3339776-03 

•3003576-37 

io 1909 6-02 

-■ 4c2 1 c4 £-0 3 

-.34P3286-05 

.b521C5t-06 

1727326-02 

-. 1977246-03 

--3402206-05 

.1235516-03 

-.1731806-02 

. 8703 756-04 

-.3112416-05 

. 1 730o7 £ -03 

16330 1 6-02 

.3357026- J3 

-.26306 Cc-05 

• 212021c--; 5 

-. 1h*5 14 E-02 

.5701 lnF-03 

-.?C9f 2CE-05 

• 2 Jc loOc-vjj 



216 

.1548346+02 

.178490E-02 

-• 88301 7E-03 

-.1194856-02 

.7033436-03 

-•151162E-05 

.2310476-05 

217 

.1555546+02 

.176529E-02 

-.3715 7CE-03 

-.915521E-03 

.7391066-03 

-•971973E-06 

.210267c-05 

216 

.1562746+02 

.1745676-02 

-• 860 122E-03 

-.6445516-03 

.6745856-03 

— • 544946E— 06 

.1731996-05 

219 

.1569946+02 

•172605E-02 

-• 848675E-03 

-.417258 E-03 

• 520503E-03 

-.2784666-06 

•Z25254C-J5 

220 

.1577156+02 

.170643E-02 

-.837227E-03 

-.262380 E-03 

.2994376-03 

-• 1 9703 7E-06 

• 730641 c-J& 

221 

. 15U436E+02 

i 1 6 666 IE— 02 

-.625 730E-03 

- • 198344 E-03 

.4249306-04 

-.2894796-06 

.235*«o6t-j6 

222 

.1591556+02 

• 166720E-02 

-.8143326-03 

-.2309906-03 

-.2150056-03 

-.56C191E-06 

-.170354E-O6 


WARNING. 

EXtRAPOLATION REQUIRED FOR » I TJ S. 



J 

TtM £» TJ(J) 

RESPONSE FUNCTION, X(, 

1 

0. 

.50666056-02 

2 

• 3000000E-01 

.39802776-02 

3 

.6000000E-01 

—•2 l 79693E-02 

4 

•90000006-01 

-• 1933 396E-01 

5 

• 12000C0E+-00 

-• 526051 5E— 01 

6 

. 15000CCE+00 

10577526*00 

7 

• 10000CCE + 00 

13040926*00 

8 

•21C0CCC6+C0 

-.27441976*00 

9 

• 24000CGE+00 

-. 383 74496*00 

10 

•27000006*00 

-.5O2O4:>0C*OO 

li 

.. 30 J0000E * 00 

-.6209243E+00 

12 

• 33000C0E ♦ 00 

— • 7 J 04 766 6 + 00 

13 

• 36000CGE + 00 

-. 0207593E*00 

14 

. 3 90 00 COE* 00 

-.03204376*00 

15 

.4200occe*oo 

- • 9055335E *00 

16 

• 45GOOCOF *00 

B822606E*00 

17 

• 48000006*00 

-. 0C883CCE *00 

18 

• 5l JOOCGE+OO 

-• 06365L 4C + 00 

19 

. 54000006*00 

-.50770296*00 

20 

. 57000006*00 

-•284 7137c *00 

21 

•6000CC06 ♦ 00 

-.20939UE-01 

22 

. 6300GCCF*00 

.2718102000 

23 

• 66 OOuOOE *00 

.58095136*00 

2* 

. 69 OOOCOt *-00 

. 69185736*00 

25 

. 72 OCOCOE *00 

.11079756*01 

26 

. 7500C'CCE*00 

. 1451 7o8 1 *0 1 

27 

. 789C0CCC *00 

. 16676346*01 

28 

• G1 JOQCOE* 00 

. 102O553E.+ O1 

29 

. 64900C0E + 00 

. 13985646*01 

30 

. 870GP COE* 00 

. 1 852 34 6 0 01 

31 

.90000006*00 

.1804 1936*01 

32 

• 93 OOOOOE *00 

. 164030 7E +01 

33 

.5600nCGE*00 

. 14094376*01 

34 

• 990Q0CCE ♦ 00 

.11231976*01 

35 

. 102 00C06* C L 

. 79461626*00 

36 

. 105C0C06+QL 

. 44 l 7595 E + 00 

37 

. 10800006*01 

. 61 C9E546-01 

38 

. Ill 00006*01 

-.27958686*00 

39 

. 1 14C00CE* 01 

-- 59040646*00 

40 

. 11700CCE*01 

— .08598926*00 

41 

. 12000 COE *01 

— . 11230896*01 

42 

. 1 23 OOCOE *01 

-. 1301C09E*01 

43 

• 1260C'C0E*0l 

-. 14135506*01 

44 

. 12900CCE*C1 

-.14575596+01 

45 

. 132 OOUOE+Ol 

-. 14524436+01 

46 

. 1350OCCE+CI 

13441076*01 

47 

. 1 3300CC6 * 01 

-.1 1909266*01 

4fl 

. l4l00C06*Cl 

L005962E +0 1 

49 

. 1440CC0E+CI 

-..7752 8786*00 

50 

. 147(JOCOE*G1 

-.521 15436*00 

51 

. 15OOOCOE*01 

-.25671 786*00 

52 

• 15 jOOCCL+OI 

.51529996-02 

5 3 

• 1 5c.00i.0t ♦ G 1 

. 2 ’>206306 *00 



54 

•15900006*01 

.47194136*00 

55 

.16200006*01 

.65637566*00 

56 

• 1650QCC6+01 

. 7984 864E*00 

57 

.16300006*01 

. 89366976*00 

58 

-17100006*01, 

• 9336542E*00 

59 

. 17400006*01 

.93364896*00 

60 

.17700006*01 

. 88237256*00 

61 

. 18000006*01 

. 7893 C89E *00 

62 

. 10300006*01 

• 6604539E+00 

63 

- 18600C0E*0i 

. 5022439E*00 

64 

. !89CuC0E*0l 

.32373456*00 

65 

. 19200006*01 

. 1337C63E *00 

66 

. 19500G0F*01 

-• 55C9901 E-Ol 

67 

. 19800006*01 

-.24593376*00 

68 

.20100006*01 

-* 41 754 io6*00 

69 

.204 00006*01 

-.56669636*00 

70 

.20700006*01 

-.6870791E+00 

71 

.21000006*01 

-.77389126*00 

72 

•21300006*01 

-.82352006*00 

73 

•21600006*01 

-. 83439996*00 

74 

•21900C0E*QL 

-. 8.0854996*00 

75 

.22200006*01 

-.74064196*00 

76 

• 22500CCE+CI 

-.65879916*00 

77 

.22800006*01 

-.54380606*00 

78 

.23100006*01 

-.41056576*00 

79 

.23400006*01 

-.26631296*00 

80 

.23700006*01 

-.11013556*00 

81 

.24000006*01 

.2699692 E -01 

82 

.24300006*01 

.16196076*00 

83 

.24600006*01 

.20123786*00 

84 

-2490OC0E*0i 

• 360265l.E*00 

85 

.25200006*01 

.45*j6940E*00 

86 

•2550UC0E+C1 

.50495886*00 

87 

.25300006*01 

• 526557 7 E* CO 

88 

.26100006*01 

.521 5 502 E *00 

39 

.26* 00006*01 

.49136152*00 

90 

.26700006*01 

.43839306*00 

91 

.27000006*01 

.36562026*00 

92 

.27300C0E*Cl 

• 2 774 53 5 E *00 

93 

• 27600 COC *01 

.17902446*00 

94 

• 2 7900C0E*01 

. 753 J 43 3E -01 

95 

.28200006*01 

-.2860634E-C1 

96 

• 28500 CCE*Cl 

-. 12770456*00 

97 

.28e00C0E*01 

-.21 755306*00 

98 

• 29 1 00C0E *01 

-.25445576*00 

99 

.29400C0E+01 

-.355i3C6E*C0 

100 

• 2970000E*0 1 

-.39775556*00 

101 

.30000006*01 

-.42055636*00 


MAX XU) - 
KIN XU) • 


18505646*01 AT TIME 
14573596*01 AT TIME 


84000006*00 

12900006*01 




* CONVOLUTION NUMBER 


3 • 




I 

FREQUENCY 

INTERPOLATED INPUT 'FREQUENCY 
RESPONSE FUNCTION 
FRF 

TRANSFORM OF 
' IMPULSE TIME 
T FRF 

THE GUST 
FUNCTION 

TRANSFORM OF THE 
TIME CONVOLUTION 
FRF • TFRF 



REAL IMAGINARY 

REAL 

IMAGINARY 

REAL IMAGINARY 


l 

0. 

-.13L924E-04 

• 333u32E-05 

-.175538E-03 

0. 

.231 577E-08 

-.5856526-09 

2 

•720158E-01 

• 462 Q79E-03 

-.373162E-02 

•074O43E-O2 

.475066E-0I 

. 1 81 3 16E-03 

— • 1C66426— 04 

3 

• 1440 32E* 00 

.500673E-03 

-.1 C1535E-01 

• 346402 6-0 1 

.8940676-01 

.9255436-03 

-.3C6915E-J3 

4 

1*2 160476 + 00 

-•611 605E-02 

-.1877816-01 

• 7505 79 E-0 1 

. 1209216*00 

.101 16CE-02 

-.2 149ol c-02 

5 

• 288063E *00 

-•1 839756-01 

-.1797616-01 

• 126 1376+00 

.1377596*00 

'• 1 5u C25E-03 

-.4802126-02 

6 

• 36007<7E*00 

-•23282 6E-01 

— • 1 10686E-0 1 

• 1829366+00 

.1372046*00 

— -.2740966-02 

-.5221 736-^2 

7 

•432O95E+0C 

-•23843 IE-01 

-.6212 12E-02 

• 240 1276*00 

.1103626*00 

- .4990056 -02 

-.4313d2t-^2 

8 

•5C4 1 10E+00 

-• 23498 6E-01 

-.34 1 599E-02 

.2919926*00 

• 6 1 53 7G6-0 1 

-.6582926-02 

-.25lJ‘ttt-02 

9 

• 576 1 26E *00 

-•2 32036E-0 1 

-.1531686-02 

.333429E*00 

.2897406-01 

-.7692366-02 

-.1103016-02 

10 

.6481 42E*0C 

-•23154 I E— C l 

15699gE-03 

.3601606*00 

-.3574806-01 

-.0344816-02 

.771 1 346-0 j 

11 

• 7201 586 >00 

-.233584E-C1 

• 106038b -02 

.3691336*00 

-• 1 G 79goC *00 

-.8507486-02 

• 29 l u32c-02 

12 

• 792 1 73E*00 

-.238845E-01 

.23424 IE-02 

.35U8 76E*00 

-.1023536*00 

— ■ 8 1 44*«3£-02 

• d 19oo 4 »6-02 

13 

.0641696*00 

-.248029E-01 

.3967096-02 

.3294 156*00 

- . 2 >34C0E +00 

-.7105176-02 

. 759l6bc-j2 

14 

• 936205E*00 

-.2614626-01 

.646273E-02 

.2324316*00 

-.3155596*00 

-.5342916-02 

• L0CG48E-ul 

15 

. 1008226*01 

-.2 7569BE-01 

.1099556-01 

.2? u992E*00 

-.3653046*00 

— • 207498E-02 

• 1 25 j jGc-jI 

16 

. 10C024E+0 1 

-.264247E-01 

• I 942 30E-0 1 

.1492706*00 

-.3985006*00 

.3795046-02 

• L 342556-4 1 

17 

.1152256*01 

-.147139E-01 

.2745046-01 

.7215456-01. 

-.4132586*00 

.1029906-01 

. 8Cc**2 1 e-02 

18 

1224276*01 

-.2 54003b- 02 

• 1 9C2 3oT -0 1 

-.51U3:>dE-02 

— . 409 1 59E *00 

.0124216-02 

.9364^26-03 

19 

• !29o2fiE*0l 

- • 1 603 09E-O2 

.1113406-01 

-.7770336-01 

-.3071026*00 

• 443 56 16 rU2 

-. 244b22 6-03 

20 

- 1360306*01 

-.3652646-02 

.6921 14E-02 

-.1411616*00 

. 34 96 3o6 *00 

• 243 54 5C-02 

• 3 CCl C0t-J3 

21 

• 1440 326* 0 L 

- • 549203E- 02 

• 5C53 746-02 

-.1919036*00 

-.249905C+00 

.257018C-02 

• 6 76d 7 7 c-Ua 

22 

. 1512336*01 

- .6 01 24 3E-02 

.417642E-02 

-.227od9£*00 

- • 2421086*00 

.25622o6-02 

. o53*lcc-v/3 

23 

-150435E*Ol 

-.7723606-02 

.374^966-02 

-.2475216*00 

-.1807026*00 

.25803oE-$2 

•46c229E-0j 

24 

• 1 65636b ♦ 0 1 

- • 63730 1l-02 

.3519416-02 

-.251 709000 

-.1200086*00 

.253C2C6-02 

• 1 L5o J4E-03 

25 

. 1 728 38 £ ♦ 0 l 

-.0841126-02 

.339768E-02 

-.2417516*00 

-.6423906-01 

•235562F-02 

-• 253**41 6-03 

26 

.1800396*01 

-.91817 It- 02 

.3331016-02 

-.2201246*00 

-. lG3913h-01 

.2075736-02 

-• 5 S2912t-03 

27 

• 187241E*0 l 

-.943202E-02 

.3 293 7 7E-02 

-• 189938 E*00 

.2118196-01 

.1722216-02 

-• 8 255o5c- J3 

28 

• 194443E * 0 l 

-.9610296-02 

.226706F-02 

154854 F*00 

.47201 1 F-0 1 

• 1 33496C-02 

— • 56 Jo 81 £— J3 

29 

• 20 1644E <-0 l 

- • 97o 73 9E-02 

.32422 IE-02 

-.1182456*00 

• 6 1 7944t-0 1 

.954551 -03 

-.9665446-03 

30 

•20084oE*0 1 

-.99056CE-02 

.3220096-02 

-.8337646-01 

• 6 5684 IE— 3 1 

.6146421-03 

-.91832 76-u3 

31 

• 2 160476*01 

— • 1 0062 9b-0l 

.3223606-02 

— .529006 F-Q l 

• 6040 &6 6- 0 1 

.3376036-03 

775-.1 76-33 

32 

• 2232 496 *01 

-•102 5346-0 1 

• 3 33 3 30E-02 

-.2871656-01 

.4U4812E-01 

.1328406-03 

-.55*0*01-03 

33 

.2304506*01 

-.1032456-01 

.3710206-02 

-.1107466-01 

.3239576-01 

.2405466-05 

— • 3 7L>*dc-0 3 

34 

■ 237652F* 01 

-.9048612-02 

• 4 1 02 3 7E-02 

-. 257132 E-02 

.1474676-01 

-.3517246-04 

-. 1557oJE-03 

35 

• 2440 546 ♦ C 1 

-.9147216-02 

.397612E-02 

-.22 75956-03 

-.210^406-02 

. 1C45326-04 

• 1 8j536C-04 

36 

« 2520 556 * 0 1 

-.8609396-02 

.3524056-02 

-.363006 6-02 

-.162-tl 76-01 

.6855616-04 

. 1270106-03 

37 

.2592576*01 

-.8101 74E-02 

.29io66E -02 

-.1116866-01 

-.2632726-01 

• 1606936-03 

« 1 82oJ5c-G3 

38 

• 2 664 5 PE*- 01 

-.7771266-02 

•21 02006-02 

-.2097b66-01 

-.3 1 678ofc- 0 1 

.2321536-03 

• 2 C Ci** uLc-u 3 

39 

• 2 73cu Gc ♦ 0 1 

-.7363726-02 

.UC771E-02 

-.3124 lot-0 1 

-.3225786-01 

• 2g 5 7o 76 -03 

. 2 02 93 1 £-03 

40 

• 28U0o 16 + 91 

— • 706 l 4ot-02 

—.4999356-03 

-.4034276-01 

-.2859696-01 

.2 7G5dCE-03 

• 222 1 0gc-c’3 

41 

• 2 RcOo 3E + 0 1 

- • 7222936-02 

-.2543876-02 

-.4703226-01 

-.2167026-01 

. 2 759 L t.E-03 

. 29*»9 79c-o3 

42 

.2^52656*01 

-.8792156-02 

-.635734E-02 

-.50524 1 E-0 1 

-. 127310E-01 

■ 3632 8CF -03 

• 43^ I3l c-<j3 

43 

.30?46gE*01 

-. 1 33530E-01 

— . 9ol 1946-02 

-.5053196-01 

-• 3 1 33o5£-02 

• 64 4 o3 OF— 03 

.52753^6-03 

44 

. 30466SF *91 

- .2 C403CE-01 

-.8 773C4E-02 

-»4724o4E-01 

• 5836626-02 

.1015176-02 

• 29>*» 1 1 c-03 

45 

• 3 16 369b + 0 1 

-.2473456-01 

-.3 852 1 56-02 

-.4125956-01 

. 131097t-01 

.1071036-02 

-. 1 63323c-j3 

46 

.324G71E*0l 

-.2547306-01 

.9690366-03 

-.3345026-01 

. 1 79366 1-0 1 

.0346976-03 

- • 4 b93 l >£-03 

47 

• 3 31 2 726 + G l 

-.2476CdE-Ul 

.4 looG56-02 

-.2484536-01 

• 14948 1F-01 

.5326341 -03 

- . 59 70u4c-u 3 


VO 

u> 



48 

•338474E+01 

- .2 39908E-01 

•6166106-02 

49 

• 3456766+01 

-.235152E-01 

.753 L63E-02 

50 

•3528776+01 

-.2342886-01 

.8645506-02 

51 

. 3600796 + 0 1 

-• 23752 5E- 0 1 

• 9 72660E-02 

52 

•3672606+01 

-.2454836-01 

• 1 094 796-01 

53 

• 3 744826+ 0 1 

- .2 59539E-01 

.1251486-01 

54 

. 3 8 1684F+ 0 1 

-.282 8076-01 

.1479846-01 

55 

.3686856+01 

- . 322406E- 01 

.1870226-01 

56 

• 3960 G7C + 01 

-.393854E-01 

•2606836-01 

57 

• 403288E* 0 1 

- .5239966-01 

• 52465 2E -0 l 

58 

« 4 L0490E + 01 

-.3709926-02 

.1376246+00 

59 

•4 17691F+01 

.7299126-01 

.9514166-02 

60 

• 424893G+01 

• 300029E -0 1 
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